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1. Introduction 

In these lecture notes we give a short introduction to the use of spectral theory in the theory 
of special functions. Conversely, special functions can be used to determine explicitly the spectral 
measures of explicit operators on a Hilbert space. The main ingredient from functional analysis that 
we are using is the spectral theorem, both for bounded and unbounded self-adjoint operators. We 
recall the main results of the theory in §2, hoping that the lecture notes become more self-contained 
in this way. For differential operators this is a very well-known subject, and one can consult e.g. 
Dunford and Schwartz 0]. 

In §3 we discuss Jacobi operators on -£^(Z>o), and we present the link to orthogonal polynomials. 
Jacobi operators on £^(Z>o) are symmetric tridiagonal matrices, and the link to orthogonal poly- 
nomials goes via the three-term recurrence relation for orthogonal polynomials. We prove Favard's 
theorem in this setting, which is more or less equivalent to the spectral decomposition of the Jacobi 
operator involved. We first discuss the bounded case. Next we discuss the unbounded case and its 
link to the classical Hamburger moment problem. This is very classical, and it can be traced back 
to at least Stone's book ^2]. This section is much inspired by Akhiezer [|l|], Berezanskii [||, §7.1], 
Deift |5|, Ch. 2] and Simon [2^], and it can be viewed as an introduction to |]l] and |^0|. Especially, 
Simon's paper [^ is recommended for further reading on the subject. See also the recent book 
|23| by Teschl on Jacobi operators and the relation to non-linear lattices, see also Deift ||5|, Ch. 2] 



for the example of the Toda lattice. We recall this material on orthogonal polynomials and Jacobi 
operators, since it is an important ingredient in §4. 

In §4 we discuss Jacobi operators on ^^(Z), and we give the link between a Jacobi operator on 
^^(Z) to two Jacobi operators on £^(Z>o) due to Masson and Repka [O], stating in particular that 
the Jacobi operator on £^(Z) is (essentially) self-adjoint if and only if the two Jacobi operators 
on £^(Z>o) are (essentially) self-adjoint. Next we discuss the example for the Meixner functions in 
detail, following Masson and Repka fl^, but the spectral measure is now completely worked out. In 
this case the spectral measure is purely discrete. If we restrict the Jacobi operator acting on ^^(Z) to 
a Jacobi operator on £^(Z>o), we obtain the Jacobi operator for the associated Meixner polynomials. 
The case of the Meixner- Pollaczek functions is briefly considered. As another example we discuss 
the second order q'-hypergeometric difference operator. In this example the spectral measure has a 



continuous part and a discrete part. Here we follow Kakehi ||l^ and [|13|, App. A]. These operators 
naturally occur in the representation theory of the Lie algebra su(l, 1), see [17|, or of the quantised 
universal enveloping algebra Uq{su{l, 1)), see ||l3|. Here the action of certain elements from the Lie 
algebra or the quantised universal enveloping algebra is tridiagonal, and one needs to obtain the 
spectral resolution. It is precisely this interpretation that leads to the limit transition discussed in 



(4.5.10) 



However, from the point of view of special functions, the Hilbert space £^(Z) is generally not 
the appropriate Hilbert space to diagonalise a second-order g-difference operator L. For the two 
examples in §4 this works nicely, as shown there. In particular this is true for the 2'/'i"Series, as 
shown in §4, see also 0] for another example. But for natural extensions of this situation to higher 



levels of special functions this Hilbert space is not good enough. We refer to |12] for the case of the 
second order g-difference operator having 3(/92-series as eigenfunctions corresponding to the big q- 
Jacobi functions, and to [15| for the case of the second order g-difference operator having gVFy-series 



as eigenfunctions corresponding to the Askey- Wilson functions. For more references and examples 
of spectral analysis of second order difference equations we refer to |14|. 



There is a huge amount of material on orthogonal polynomials, and there is a great number of 
good introductions to orthogonal polynomials, the moment problem, and the functional analysis 
used here. For orthogonal polynomials I have used P], P, Ch. 2] and |p4|. For the moment problem 



there are the classics by Shohat and Tamarkin |19] and Stone [^2|, see also Akhiezer ||l|], Simon 
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|20| and, of course, Stieltjes's original paper [21| that triggered the whole subject. The spectral 
theorem can be found in many places, e.g. Dunford and Schwartz [^ and Rudin [18|. 

Naturally, there are many more instances of the use of functional analytic results that can be 
applied to special functions. As an example of a more qualitative question you can wonder how 
perturbation of the coefficients in the three-term recurrence relation for orthogonal polynomials 
affects the orthogonality measure, i.e. the spectral measure of the associated Jacobi operator. 
Results of this kind can be obtained by using perturbation results from functional analysis, see e.g. 
Dombrowski [^ and further references given there. 

Acknowledgement. I thank the organisors, especially Renato Alvarez-Nodarse and Francisco Mar- 
cellan, of the summer school for inviting me to give lectures on this subject. Moreover, I thank the 
participants, as well as Johan Kustermans, Hjalmar Rosengren and especially Wolter Groenevelt, 
for bringing errors to my attention. 
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2. The spectral theorem 

2.1. Hilbert spaces and bounded operators. 

(2.1.1) A vector space Ti. over C is an inner product space if there exists a mapping (•, •) : TCxTC -^ C 
such that for all u,v,w € TC and for all a,b £ C we have (i) {av + bw,u) = a{v,u) + b{w,u), (ii) 
{u,v) = {v,u), and (iii) {v,v) > and {v,v) = if and only if v = 0. With the inner product 
we associate the norm \\v\\ = ||w||?^ = y^{v,v), and the topology from the corresponding metric 
d{u,v) = \\u — v\\. The standard inequality is the Cauchy-Schwarz inequality; |(ii, v)| < ||u||||f ||. A 
Hilbert space W is a complete inner product space, i.e. for any Cauchy sequence {x„}„ in 7i, i.e. 
Ve > 3N G N such that for all n,m > N ||x„, — Xm\\ < £> there exists an element x £ 7i such that 
Xn converges to x. In these notes all Hilbert spaces are separable, i.e. there exists a denumerable 
set of basis vectors. 

(2.1.2) Example. ^^(Z), the space of square summable sequences {ak}keZi and £^(Z>o), the 
space of square summable sequences {ak}kez>o^ ^^'^ Hilbert spaces. The inner product is given by 
({flfcji {bk}) = X^fc '^kbk- An orthonormal basis is given by the sequences e^ defined by {ek)i = 6k,i, 
so we identify {ofc} with ^f^akCk- 

(2.1.3) Example. We consider a positive Borel measure fi on the real line M such that all moments 
exist, i.e. J^ \x\"^ dfi{x) < oo for all m G Z>o. Without loss of generality we assume that // is a 
probability measure, f^ dn{x) = 1. By L?'{n) we denote the space of square integrable functions on 
M, i.e. Jjg |/(x)p dfj,{x) < oo. Then L'^in) is a Hilbert space (after identifying two functions / and g 
for which f^ \f{x) — g{x)\'^ dfi{x) = 0) with respect to the inner product {f,g) = f^ f{x)g{x) dix{x). 
In case /x is a finite sum of discrete Dirac measures, we find that L'^{h) is finite dimensional. 

(2.1.4) An operator T from a Hilbert space 7i into another Hilbert space K, is linear if for all 
u,v £ TC and for all a, 6 G C we have T{au + bv) = aT{u) + bT{v). An operator T is bounded if 
there exists a constant M such that ||Tn||yc < M||m||7^ for all u £ TC. The smallest M for which this 
holds is the norm, denoted by ||T||, of T. A bounded linear operator is continuous. The adjoint 
of a bounded linear operator T : TC ^ K, \s a. map T* : K, ^ TC with {Tu,v)ic = {u,T*v)'n. We 
call T: TC ^ TC self-adjoint if T* = T. It is unitary if T*T = 1->y and TT* = 1a:. A projection 
P : TC ^ TC IS a, linear map such that P^ = P. 

(2.1.5) We are also interested in unbounded linear operators. In that case we denote {T,'D(T)), 
where T){T), the domain of T, is a linear subspace of TC and T: T>[T) -^ TC. Then T is densely 
defined if the closure of 2?(T) equals TC. All unbounded operators that we consider in these notes 
are densely defined. If the operator (T — z), z G C, has an inverse R{z) = {T — z)~^ which is densely 
defined and is bounded, so that R{z), the resolvent operator, extends to a bounded linear operator 
on TC, then we call z a regular value. The set of all regular values is the resolvent set p{T). The 
complement of the resolvent set p{T) in C is the spectrum cr{T) of T. The point spectrum is the 
subset of the spectrum for which T — z is not one-to-one. In this case there exists a vector v £ TC 
such that (T — z)v = 0, and z is an eigenvalue. The continuous spectrum consists of the points 
z G cr{T) for which T — z is one-to-one, but for which (T — z)TC is dense in TC, but not equal to 
TC. The remaining part of the spectrum is the residual spectrum. For self-adjoint operators, both 



bounded and unbounded, see |(2.3.3)| , the spectrum only consists of the discrete and continuous 
spectrum. 

(2.1.6) For a bounded operator T the spectrum cr{T) is a compact subset of the disk of radius 
||r||. Moreover, if T is self-adjoint, then cr{T) C M, so that cr{T) C [— ||r||, ||r||] and the spectrum 
consists of the point spectrum and the continuous spectrum. 
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2.2. The spectral theorem for bounded self-adjoint operators. 

(2.2.1) A resolution of the identity, say E, of a Hilbert space TC is a projection valued Borel 
measure on M such that for all Borel sets ^, i? C M we have (i) E{A) is a self-adjoint projection, (ii) 
E{AnB) = E{A)E{B), (iii) E(0) = 0, E{R) = In, (iv) AnB = implies E{AUB) = E{A)+E{B), 
and (v) for all u,v £ 7i the map A i-^ E^^vi-^) = {E{A)u,v) is a complex Borel measure. 

(2.2.2) A generalisation of the spectral theorem for matrices is the following theorem for bounded 
self-adjoint operators, see @, §X.2], |l|, §12.22]. 

Theorem. (Spectral theorem) Let T: TC ^ 7i be a bounded self-adjoint linear map, then there 
exists a unique resolution of the identity such that T = J^tdE{t), i.e. {Tu,v) = J^tdEu^vit)- 
Moreover, E is supported on the spectrum (t{T), which is contained in the interval [— ||T||, ||T||]. 
Moreover, any of the spectral projections E{A), A <ZM. a Borel set, commutes with T . 

A more general theorem of this kind holds for normal operators, i.e. for those operators satisfying 

(2.2.3) Using the spectral theorem we define for any continuous function / on the spectrum 
a{T) the operator f{T) by /(T) = j^f{t)dE{t), i.e. {f{T)u,v) = J^fit)dE^,,{t). Then /(T) 
is bounded operator with norm equal to the supremum norm of / on the spectrum of T, i.e. 
||/(T)|| = sup^.go.(2-) 1/(3^)1- This is known as the functional calculus for self-adjoint operators. In 
particular, for z G p{T) we see that /: x i— > (x — z)~^ is continuous on the spectrum, and the 
corresponding operator is just the resolvent operator R{z) as in (2.1.5)| . The functional calculus 



can be extended to measurable functions, but then ||/(r)|| < su])^^^r']^\ |/(a;)|. 

(2.2.4) The spectral measure can be obtained from the resolvent operators by the Stieltjes-Perron 
inversion formula, see [0, Thm. X.6.1]. 

Theorem. The spectral measure of the open interval (a, b) cM is given by 



I rb-S 



Eu y( (a, b)) =liinliin / (R(x + ie)u,v) — (R(x — ie)u,v) dx. 

' 6io £iO 27^^ J^+s 

The limit holds in the strong operator topology, i.e. T„x — > Tx for all x £ 7i. 

2.3. Unbounded self-adjoint operators. 

(2.3.1) Let {T,V{T)), with P(r) the domain of T, be a densely defined unbounded operator on 
7i, see K2.1.5)| . We can now define the adjoint operator (T*,2?(T*)) as follows. First define 

V{T*) = {u S H I n 1-^ {Tu,v) is continuous on P(T)}. 

By the density of P(T) the map u i— > {Tu, v) for v G T>[T*) extends to a continuous linear functional 
uj: 7i — > C, and by the Riesz representation theorem there exists a unique w G H such that 
uj{u) = (n, w) for all u G TC. Now the adjoint T* is defined by T*v = w, so that 

{Tu,v) = {u,T*v) yueV{T),yveV{T*). 

(2.3.2) If T and S are unbounded operators on TC, then T extends S, notation S C T, if 
X'(S') C 'D(T) and Sv = Tv for all v £ T>{S). Two unbounded operators S and T are equal, S = T, 
if S C T and T C S, or S and T have the same domain and act in the same way. In terms of the 
graph 

g{T) = {{u,Tu) I u G V{T)} cTCxTC 
we see that 5 C T if and only if G{S) C GiT). An operator T is closed if its graph is closed in the 
product topology of TC x TC. The adjoint of a densely defined operator is a closed operator, since 
the graph of the adjoint is given as 

giT*) = {{-u,Tu)\uGV{T)}^, 
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for the inner product ((n, v), {x, y)) = {u, x) + {v, y) on TC x 7i, see [18, 13 



(2.3.3) A densely defined operator is symmetric if T C T*, or, using the definition in (2.3.1) 

{Tu,v) = {u,Tv), y u,v ^V{T). 

A densely defined operator is self-adjoint if T = T*, so that a self-adjoint operator is closed. The 
spectrum of an unbounded self-adjoint operator is contained in M. Note that V{T) C P(T*), so 
that T>{T*) is a dense subspace and taking the adjoint once more gives {x** ^T>{T**)) as the minimal 
closed extension of {T,T){T)), i.e. any densely defined symmetric operator has a closed extension. 
We have T C T** C T* . We say that the densely defined symmetric operator is essentially self- 
adjoint if its closure is self-adjoint, i.e. if T C T** = T*. 

(2.3.4) In general, a densely defined symmetric operator T might not have self-adjoint extensions. 
This can be measured by the deficiency indices. Define for z G C\M the eigenspace 

N^ = {v £V{T*) I T*v = zv}. 

Then dimA'^j; is constant for 9z > and for Qz < 0, 0, Thm. XII. 4. 19], and we put n_|- = dim Aj 
and n_ = dimA^_j. The pair (n^,n-) are the deficiency indices for the densely defined symmetric 
operator T. Note that if T commutes with complex conjugation of the Hilbert space TC then we 
automatically have n_|_ = n_. Note furthermore that if T is self-adjoint then n+ = n_ = 0, since a 
self-adjoint operator cannot have non-real eigenvalues. Now the following holds, see |Q, §XII.4]. 

Proposition. Let {T,T){T)) he a densely defined symmetric operator. 

(i) T){T*) = T>{T**) (B Ni (B N-i, as an orthogonal direct sum with respect to the graph norm for 

T* from {u,v)t* = {u,v) + {T*u,T*v). As a direct sum, T){T*) = T){T**) + N^ + N^ for general 

z £ C\M. 

(ii) Let U be an isometric bijection U : Ni ^ N^i and define {SjD^S)) by 

V{S) = {u + v + Uv\ue V{T**), V G Ni}, Sw = T*w, 

then {S,'D{S)) is a self-adjoint extension of {T,V{T)), and every self-adjoint extension ofT arises 
in this way. 

In particular, T has self-adjoint extensions if and only if the deficiency indices are equal; n+ = n_. 
T** is a closed symmetric extension of T. We can also characterise the domains of the self-adjoint 
extensions of T using the sesquilinear form 

B{u, v) = {T*u, v) - {u, T*v), u,v£ V{T*), 

then V{S) = {u e V{T*) \ B{u,v) = 0, Vw G ^{3)}. 

2.4. The spectral theorem for unbounded self-adjoint operators. 



(2.4.1) With all the preparations of the previous subsection the Spectral Theorem (2.2.2) goes 
through in the unbounded setting, see [0, §XII.4], |l^, Ch. 13]. 

Theorem. (Spectral theorem) Let T: T){T) -^ TC be an unbounded self-adjoint linear map, then 
there exists a unique resolution of the identity such that T = J^tdE{t), i.e. {Tu,v) = J^tdEu^vit) 
for u G D^T), V £ TC. Moreover, E is supported on the spectrum cr(T), which is contained in M. 
For any bounded operator S that satisfies ST C TS we have E{A)S = SE{A), A C M a Borel set. 



Moreover, the Stieltjes- Perron inversion formula (2.2.4) remains valid; 



1 f^^^ 

EuvdcL^bX) =\mi\mi / (Rix + ie)u,v) — {R(x — ie)u,v) dx. 

' ^' '^ Sio ao 2-711 J^+s 
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(2.4.2) As in (2.2.3)| we can now define f{T) for any measurable function / by 



{f{T)u,v)= / f{t)dEuAt), ueV{f{T)), ven, 

Jm 

where V{f{T)) = {u e H \ j^ |/(t)P dEu,u{t) < 00} is the domain of /(T). This makes /(T) into 
a densely defined closed operator. In particular, if / G L°°{U.), then f{T) is a continuous operator, 
by the closed graph theorem. This in particular applies to f{x) = {x — z)~^ , z G /o(T), which gives 
the resolvent operator. 



8 erik koelink 

3. Orthogonal polynomials and Jacobi operators 

3.1. Orthogonal polynomials. 

(3.1.1) Consider the Hilbert space L'^{^) as in Example |(2.1.3)| . Assume that ah moments exist, 
so that all polynomials are integrable. In applying the Gram-Schmidt orthogonalisation process to 
the sequence {1, x, x'^,x^, . . . } we may end up in one of the following situations: (a) the polynomials 
are linearly dependent in L^(/i), or (b) the polynomials are linearly independent in L^(/i). In case 
(a) it follows that there is a non-zero polynomial p such that J^\p{x)\'^ dfi{x) = 0. This implies 
that /i is a finite sum of Dirac measures at the zeros of p. From now on we exclude this case, but 
the reader may consider this case him/herself. In case (b) we end up with a set of orthonormal 
polynomials as in the following definition. 

Definition. A sequence of polynomials {pn}^=o with deg(p„) = n is a set of orthonormal polyno- 
mials with respect to fi if J^pn{x)pm{x) dfj,{x) = 5n,m- 

Note that the polynomials pn are real- valued for x E M, so that its coefficients are real. Moreover, 
from the Gram-Schmidt process it follows that the leading coefficient is positive. 

(3.1.2) Note that only the moments nik = j^x^ d^{x) of // play a role in the orthogonalisation 
process. The Stieltjes transform of the measure fi defined by w{z) = J^{x — z)~^ dfi{x), z G C\M, 
can be considered as a generating function for the moments of [i. Indeed, formally 

z ]^\-xlz z ^^^J^ z ^^^z 

In case supp(^) C [—A, A] we see that \mk\ < 2A'' implying that the series in (|3.lD is absolutely 
convergent for |2:| > A. In this case we see that the Stieltjes transform w{z) of fj, is completely 
determined by the moments of /i. In general, this expansion has to be interpreted as an asymptotic 
expansion of the Stieltjes transform w{z) as \z\ -^ oo. We now give a proof of the Stieltjes inversion 



formula, cf. |(2.2.4) 



Proposition. Let fi be a probability measure with finite moments, and let uj{z) = J^{x — z) ^ diJ,{x) 
be its Stieltjes transform, then 

1 /''' 11 

lim - / Q{w{x + ie)) dx = fi{{a, b)) + -^({a}) + -fJ-i{b}). 

£lO IT J a 2 2 

Proof. Observe that 



2i'^(w(z)) = w(z) — w{z) = w{z) — w{z) = / z dfi(x) 

Jrx - z X - z 

z — z f -^ z 
77jdfi(x) = 2i / 77^dfi(x), 

-^..x-z\^ J^\x-z\^ 

so that 

^^"^" + ''^^ = L \s-{x + ^e)\- ^^^^) = L {s-xY + e- ^^^^^^ 
Integrating this expression and interchanging integration, which is allowed since the integrand is 
positive, gives 

'^{wix + ie\) dx = I I -7^ T^dxduis). (3.2) 

The inner integration can be carried out easily; 



b g r{b-s)/e ]_ 



Xe{s) = / 7 Tg- — 2^^= / T"^ — ^dy = arctany 

{s-xy+£^ J(^a-s)/e 1 + r 



{b-s)/e 
(a—s)/e 
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hy y = {x — s)/e. It follows that < Xe(s) < vr, and 

,. , . I vr, for a < s < 6, 

e[0 I 2vr, tor s = o or s = 0. 

It suffices to show that we can interchange integration and the limit e | in ( |3.2D . This follows from 
Lebesgue's dominated convergence theorem since /i is a probability measure and < Xe{s) < vr. D 



As a corollary to the proof we get, cf. (2.2.4) , (2.4.1) , 

rb-5 



1 f 
limlim— / '^(wix + ie)] dx = n((a,h)] . 



40 5iO vr 7„+5 
We need the following extension of this inversion formula in (3.3.10)| . For a polynomial p with 



real coefficients we have 



If f 11 

lim— / Q(p{x + ie)w{x + ie)) dx = / p{x) dfi{x) + -p{a) fi{{a}) + -p{b) fi{{b}) . (3.3) 

£■10 TT Ja J{a,b) 2 2 

We indicate how the proof of the proposition can be extended to obtain (|3.3|) . Start with 



r {s-x)^{p{x + ie))+e^{p{x-ie)) 
''-j[p{x + ie)w{x + ie)) = / ^ — -. -^ tt" -dfj.{s). 

Jr [s-xy + e^ 

Integrate this expression with respect to x and interchange summations, which is justified since 
^(^(x + ie)) and 3f?(p(x — is)) are bounded on (a, b). This time we have to evaluate two integrals. 
The first integral 

^ {s-x)Q{p{x + ie)) _ /■C'-'')/^ -y9(p(ey + s + ie)) 

(s-x)2+e2 ^ J(a-s)/e 2/^ + 1 

can be estimated, using 9(p(x + ie)) = 0{e) uniformly on [a, b], by 

eM / -^Ji^dy = eM\nyJy^ + l 

J{a-s)/e y + 1 {a~s)/e 

This term tends to zero independently of a, 6 and s, since eln(-2- + 1) = — 2elne + £\n{A^ + £^) 
which can be estimated by O(elne) with a constant independent of A. The other integral can be 
dealt with as in the proof of the proposition. Next Lebesgue's dominated convergence theorem can 
be applied and ( |3.3| ) follows. 

(3.1.3) The following theorem describes the fundamental property of orthogonal polynomials in 
these notes. 

Theorem. (Three term recurrence relation) Let {pk\^=Q be a set of orthonormal polynomials in 
L^{fj,), then there exist sequences {ofclfcLo' i^fclfcLo' with a^ > and 6^ G M, such that 

xpk{x) = akPk+i{x) + bkPk{x) + ak-iPk~i{x), k>l, (3.4) 

xpq{x) = aQpi{x)+bQPQ{x). (3.5) 

Moreover, if fi is compactly supported, then the coefficients Ok and bk are bounded. 



Note that (3.4), ( |3.5| ) together with the initial condition po{x) = 1 completely determine the 
polynomials Pk{x) for all A; G N. 

Proof. The degree of xpk{x) is A; + 1, so there exist constants Cj such that xpk{x) = Ylii=Q CiPi{x). 
By the orthonormality properties of pk it follows that 

Pi{x)xpk{x)dfi{x). 
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Since the degree of xpi{x) is i + 1, we see that Cj = for i + 1 < fc. Then 

,2 



bk = Ck= / x{pk{x))) dfi{x) G M. 

Moreover, c^+i = J^Pk+i[x)xpk{x) d^{x) and Cfc_i = f^pk-i{x)xpk{x) dfi{x) display the required 
structure for the other coefficients. The positivity of a^ fohows by considering the leading coefficient. 
For the last statement we observe that 



Iflfcl = I / xpk+i{x)pk{x) dfi{x)\ < / \pk+iix)\\pkix)\ dfi{x) sup |x| 

Jr Jr x&s\ipp{fi) 

<\\Pk+l\\mp)\\Pk\\L2(^,) sup |x| = sup |x| < OO, 

x£s\ipp{fi) x&s\ipp{fi) 

since ||pfc||L2(/t) = 1 and supp(/u) is compact. In the second inequality we have used the Cauchy- 
Schwarz inequality (2.1.f)| . Similarly, 

\bk\ < l|PA:|li2(^) sup |x| = sup |x| < OO 
x(^snpp{fi) a;£supp(/i) 

gives the estimate on the coefficients bk- □ 



(3.1.4) We observed that (|3.4| ) and ( ^.51) together with an initial condition for the degree zero 
component completely determines a solution for the recurrence (|3.4| ), ( |3.5D . We can also generate 
solutions of (|3.4D by specifying the initial values for k = and k = 1. Prom now on we let ^^(a;) 
be the sequence of polynomials generated by (|3.4| ) subject to the initial conditions ro{x) = and 
ri(x) = Og" . Then r^ is a polynomial of degree /c — 1, and (3.5) is not valid. The polynomials 
{^k}'kLo ^^^ t^6 associated polynomials. 

Lemma. The associated polynomial rj. is given by 

, . f Pk{x)-Pk{y) ^ , . 

rk{x) = / did{y). 

Jr x-y 

Proof. It suffices to show that the right hand side, denoted temporarily by Qkix), satisfies the 
recurrence (3.4) together with the initial conditions. Using Theorem (3.1.3) for Pk{x) and the 
definition of qkix) we obtain 

xqk{x) =akqk+i{x) + bkQkix) + ak-iQk-i 



+ 



f akPk+i{y) + bkPkjy) + ak-iPk-i{y) - xpkjy) , . . 
Jr x-y 



Using Theorem (3.1.3) again shows that the integral equals — jmPk{y) dfj.{y), which is zero for k > 1 
and —1 for A; = by the orthogonality properties. Hence, ( p.4p is satisfied. Using po{x) = 1 we 
find qo{x) = and using pi{x) = a^ (x — 6o) gives qi{x) = a^ . D 



Considering 



X- Z Jk X 

immediately proves the following corollary. 



Pkix) , , . f Pk{x)-Pk{z) , . . . s , - . t \ 

dfi{x) = / dfi(x)+pk{z) / dfi{x) 



Corollary. Let z G C\M be fixed. The k-th coefficient with respect to the orthonormal set {pk}^=Q 
in L'^ifJ.) of x ^^ {x — z)^^ is given by w{z)pk{z) + rk{z). Hence, 

OO „ 

'S2\w{z)pk{z) +rk{z)\'^ < / \x - z\^'^ dij.{x) < OO. 

T^r. JR. 



k=0 
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The inequality follows from the Bessel inequality. If the {pn}^=o is an orthonormal basis of L^(/x) 
then we can write an equality by Parseval's identity. 

(3.1.5) Since rfc(y) is another solution to (^^), multiplying ( |3.4D by ri:{y) and ( |3.4| ) for rjt(y) by 
Pk{x), subtracting leads to 

{x - y)Pk{x)rk{y) =ak{pk+i{x)rk{y) - Pk{x)rk+i{y)) 

(3-6) 
- ak-i{pk{x)rk-i{y) - Pk-i{x)rk{y)) 

for A; > 1. Taking x = y vn ( |3.6D , we see that the Wronskian, or Casorati determinant, 

[p,r\k{x) = ak{pk+i{x)rk{x) - Pk{x)rk+i{x)) (3.7) 

is independent of A; € Z>o, and taking A; = gives [p,r\k[x) = [p, r] = —1. This also shows that pk 
and Tfc are linearly independent solutions to ( |3.4] ). 

On the other hand, replacing r^ by pk and summing we get the Christoffel-Darboux formula 

n-l 

{X - y)^Pk{x)pk{y) = an-l{pn{x)Pn-l{y) - Pn-l{x)Pn{y)) ■ (3-8) 

fc=0 

The case x = y \s obtained after dividing hy x — y and in the right hand side letting y ^ x. This 
gives 

n-l 

"^Pkixf = an-l{Pn{x)Pn-l{x) - p'n~l{x)Pn{x)) ■ 

fc=0 

3.2. Moment problems. 

(3.2.1) The moment problem consists of the following two questions: 

1. Given a sequence {mQ,mi,m2, ■ ■ ■}, does there exist a positive Borel measure /i on R such 
that rrafc = J x^ dfi{x)7 

2. If the answer to problem 1 is yes, is the measure obtained unique? 

Note that we can assume without loss of generality that tuq = 1. This is always assumed. 

(3.2.2) In case supp(;u) is required to be contained in a finite interval we speak of the Haussdorf 
moment problem (1920). In case supp(^) is to be contained in [0, oo) we speak of the Stieltjes 
moment problem (1894). Finally, without a condition on the support, we speak of the Hamburger 
moment problem (1922). Here, a moment problem is always a Hamburger moment problem. The 
answer to question 1 can be given completely in terms of positivity requirements of matrices com- 



posed of the mj's, see Akhiezer [y, Shohat and Tamarkin ||T^, Simon |2y], Stieltjes |^, Stone ||22 |. 
In these notes we only discuss an answer to question 2, see §3.4. 

In case the answer to question 2 is affirmative, we speak of a determinate moment problem and 
otherwise of an indeterminate moment problem. So for an indeterminate moment problem we have 
a convex set of probability measures on M solving the same moment problem. The fact that this 



may happen has been observed first by Stieltjes |21]. The Haussdorf moment problem is always 



determinate as follows from (3.1.2) and the Stieltjes- Perron inversion formula, see Proposition 



WM. 



For a nice overview of the early history of the moment problem, see Kjeldsen |11|. 
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3.3. Jacobi operators. 
(3.3.1) A tridiagonal matrix of the form 



J 



/bo ao ...\ 

aQ bi ai . . 

ai 62 02 .. 

a2 &3 03 . . 

V; ••. ••. ••. ••. J 

is a Jacobi operator, or an infinite Jacobi matrix, if 6j G M and Cj > 0. If ai = for some i, the 
Jacobi matrix sphts as the direct sum of two Jacobi matrices, of which the first is an (i + 1) x (i + 1)- 
matrix. 

We consider J as an operator defined on the Hilbert space ^^(Z>o), see Example [2.1.2) . So with 



respect to the standard orthonormal basis {ek}kez>o of £ (Z>o) the Jacobi operator is defined as 



J Gk = < 

[ooei + boeo, 



k>l, 
k = 0. 



(3.' 



Note the similarity with Theorem (3.1.3) . So to each probability measure on R with finite moments 
we associate a Jacobi operator on the Hilbert space £^(Z>o) from the three-term recurrence relation 
for the corresponding orthonormal polynomials. However, some care is necessary, since ( |3.9[ ) might 
not define a bounded operator on ^^(Z>o). 

From (3.9) we extend J to an operator defined on 2?(Z>o), the set of finite linear combinations 
of the elements Ck of the orthonormal basis of £^(Z>o). The linear subspace P(Z>o) is dense in 
£2(Z>o). From (El) it follows that 



{Jv,w) = {v,Jw), y V,W £ ViZyo), 



(3.10) 



so that J is a densely defined symmetric operator, see (2.3.3) . In particular, if J is bounded on 

P(Z>o), J extends to a bounded self-adjoint operator by continuity. 

Lemma. (3.3.2) e^ = Pk{J)eo for some polynomial Pk of degree k with real coefficients. In 

particular, cq is a cyclic vector for the action of J, i.e. the linear subspace {J cq \ k G Z>o} is 

dense in £^(Z>o). 

Proof. It suffices to show that e^ = Pk{J)eQ for some polynomial of degree k, which follows easily 
from (|3.9|) and Oj > using induction on k. D 



Lemma. (3.3.3) If the sequences {ok} and {bk} are bounded, say sup^ jafcl + sup^ \bk\ < M < 00, 
then J extends to a bounded self-adjoint operator with \\J\\ < 2M . On the other hand, if J is 
bounded, then the sequences {ok} and {bk} are bounded. 



Proof. If {ok}, {bk} are bounded, then, with v = YlT=o'^f^^k ^ -^(^ 



'>o; 



1, 



I'^'^ll^ = X^ \akVk-i + bkVk + flfc-i^fc+iP- 



fc=0 



Let sup^ \ak\ = A, sup^r. \bk\ = B with A + B < M, then each summand can be written as 

al\vk-i\'^ + bllvkl"^ + al^ilvk+il"^ + 2akbk^{vk-iVk) + 2akak-i?fi{vk-iVk+i) + 2bkak^{vk+iVk) 
<A\\vk-i\^ + hfc+il') + B^\vk\^ + 2A^mvk-iWTl)\ + 2AB\^{vk-ivj:)\ + 2ABMvk+iv]:)\. 
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Using the bounded shift operator 5: e^ i-^ e^+i we have, using the Cauchy-Schwarz inequahty 



(2.1.1)1 and ||5|| = 1, 

\\Jvf <2A'^ + B^ + 2A'^\{S^v,v)\+AAB\{Sv,v)\ < AA^ + B^ + AAB 

= {A + Bf + 2AB + 3^2 = 2{A + Bf + 2A^ - B"^ < AM^. 

By continuity J extends to a bounded operator on £^(Z>o). 

To prove the reverse statement, we have |(JeA:,e/)| < ||J|| implying that \ak\ < ||J|| (take 
/ = A; + 1) and \bk\ < \\J\\ (take I = k). D 

(3.3.4) Assume J is bounded, then J is self-adjoint and we can apply the spectral theorem for 



bounded self-adjoint operators, see Theorem (2.2.2) . Thus 



{Jv,w)= / tdE^^^{t), yv,w£i'^{Z>o). 

In particular, we define the measure /x(A) = £'eo,eo(^) = {E{A)eo,eo). Since £■ is a resolution of 
the identity, E{A) is an orthogonal projection implying that // is a positive Borel measure. Indeed, 
H{A) = {E{Afeo,eo) = {E{A)eo,E{A)eQ) > for any Borel set A. Moreover, ^(M) = /, so that 
^ is a probability measure. The support of /x is contained in the interval [— 1|^||, ||^||], since J is 
a bounded self-adjoint operator. In particular, fi has finite moments. Hence the spectral theorem 
associates to a bounded Jacobi operator J a compactly supported probability measure fj,. Moreover, 
the spectral measure E is completely determined by fi. Indeed, 

{E{A)ek,ei) = {E{A)Pk{J)eo,Pi{J)eo) 

= {PiiJ)Pk{J)E{A)eo,eo)= [ Pk{x)Pi{x)df,{x), ^^'^^^ 

Ja 



where the polynomials Pk are as in Lemma (3.3.2) using the self-adjointness of J and the fact that 

the spectral projections commute wit^ ,^y , 

Theorem. (3.3.5) Let J be a bounded Jacobi operator, then there exists a unique compactly sup- 
ported probability measure fi such that for any polynomial P the map U : P(J)eo >-^ P extends to a 
unitary operator ^■^(Z>o) -^ -^^(/^) with UJ = MU, where M: Lp'{^) -^ L'^ilj) is the multiplication 
operator [M f){x) = xf{x). Moreover, let p^ = Ue^, then the set {pk}'k'=o ^■^ ^^^ ^^^ ^f orthonormal 
polynomials with respect to n; 

Pk{x)pi{x) dfi{x) = 5k,i. 



Proof. By Lemma (3.3.2) we see that U maps a dense subspace of i {Z>q) onto a dense subspace 
of i^(/i), since the polynomials are dense in Lp'ip) because /i is compactly supported. Using ( 3.11 ) 
we see that for any two polynomials P, Q we have 

(P(J)eo,Q(J)eo) = (Q(J)P(J)eo,eo) = / Q{x)P{x)d^{x) 

JR 

= {P,Q)lH^) = (t/P(J)eo,C/g(J)eo)L2(^), 

or U is unitary, and it extends uniquely to a unitary operator. 

To show that U intertwines the Jacobi operator J with the multiplication operator we show that 
UJck = MU Cfc for all k G Z>o. Define p^ = Ue^, then p^ E L'^ilJ-) is a polynomial of degree k and 
the set {pk}'kLo ^^ ^^^ ^^^ °^ orthonormal polynomials with respect to fi; 

Pk{x)pi{x) dfi{x) = (C/efc,t/ei)i2(^) = (efc,e/) = 5k,i 
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by the unitarity of U. It now suffices to show that the coefficients in the three-term recurrence 



relation for the orthogonal polynomials as in Theorem (3.1.3) correspond to the matrix coefficients 



of J. This is immediate using the functional calculus of (2.2.3) , (3.11) and the explicit expressions 
for the coefficients in the three-term recurrence relation given in the proof of Theorem (3.1.3) ; 



ak = {Jek,ek+i) = / x dE^^ ^ek+ii^) = / xpkix)pk+iix) dfi{x), 

JR JR 

bk = {Jek,ek) = / X dEe^^^Ckix) = / x{pk{x)) dfi{x). 

JR JR 

To show uniqueness, we observe that the moments of /x are uniquely determined by the fact that 
{Pk}kL(} is a set of orthonormal polynomials for L?'{^). Since the measure is compactly supported 
its Stieltjes transform is analytic in a neighbourhood of oo. Using the Stieltjes inversion formula 



of Proposition |(3.1.2)| , we see that the compactly supported measure is uniquely determined by its 
moments. D 



Theorem (3.3.5) is called Favard's theorem restricted to the case of bounded coefficients in the 
three-term recurrence operator. It states that any set {pfcjfcLo "-"^ polynomials generated by ( p.4p , 
( |3.5| ) with the initial condition Po{x) = 1 with o^ > 0, 5^ G M are orthonormal polynomials with 
respect to some, see §3.2, positive probability measure /x. 

Corollary. The moment generating function for fj, is in terms of the resolvent R{z) for the Jacohi 
operator J ; 

= {R{z) eo,eo) = {{J - z) eo,eo), ze 



(3.3.6) The asymptotically free solution to J f{z) = zf{z) for z G C\M, is the element 
f{z) = {fk{z)}'^^Q satisfying (^Jf{z))^ = z fk{z) for A; > 1 (and in general not for A: = 0) and 
Z^fcLo \fk{z)\'^ < cxD. The asymptotically free solution encodes the spectral measure of J, hence the 
probability measure /U of Theorem (3.3.5) . 

Proposition. Let J be a bounded Jacobi matrix. Take z £ C\M fixed. Then f{z) = (J — z)~^eQ is 
the asymptotically free solution for the Jacobi operator J. There exists a unique w{z) £ C\M such 
that fk{z) = w{z)pk{z) +rk{z), with pk, rk the polynomials as in [3.1.3), (3.1.4)| . Moreover, w is 
the Stieltjes transform of the measure /u. 



Proof. We have already observed in Corollary (3.1.4)| that {w{z)pk{z) +rk{z)}'^^Q G ^^(Z>o). Next 
we consider 



fk{z) = {{J - z) eo,efc) 



Pk{x} 



dfi{x) = w{z)pk{z) + rk{z) 



___ x — z 
again by |(3.1.4)| . Hence, {Jf{z))^ = z fk{z) for k>l. 

It remains to show uniqueness. If not, then there would be two linearly independent solutions, 
so that we could combine to get YlT=o \Pk{z)\^ < oo for z G C\R, so that J would have a non-real 
eigenvalue z contradicting the self-adjointness of J. D 



(3.3.7) Note that we need linix- 
This implies 



w{z) 



,w{z)Pk{z)+rk{z) 

r rk{z) 
= - hm — --, 

k^oopk{z) 



in order to have {A(z)}^o ^ ^^(^ 



->o) 



z G 



assuming that the limit in the right hand side exists. The fraction 
to the following lemma. 

Lemma. The zeroes of pk(x) are real and simple. 



Pk{z) 



has no non-real poles due 
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Proof. Define 



(bo 


«o 








\ 


flo 


h 


ai 















ajv-2 


bN-1 aN~i 


^0 









aN-i bN J 



Jn 



i.e. a tridiagonal matrix that is obtained from J by keeping only the first (A^ + 1) x (A'^ + 1) block 
matrix. Then J^r = J^, and it follows that its spectrum is real. Moreover, its spectrum is simple. 
Indeed, ( Jat — A)/ = and /o = implies /i = and hence fk = 0. So if the multiplicity of the 
eigenspace is more than one, we could construct a non-zero eigenvector with /o = 0, a contradiction. 
On the other hand we have, from Theorem (3.3.5), 



{Jn - z) 



so that the eigenvalues of Jn are the zeroes oipN+i- 



/po{z)\ 

Pi{z) 


= -aN 


( ° ^ 



\PNiz)J 




\PN+l{z)J 



D 



(3.3.8) The norm of the asymptotically free solution for fixed z 6 C\M can be expressed in terms 
of the Stieltjes transform of /i; 



^\w{z)pk{z) +rk{z)\^ = {{J - z) ^eo,{J-z) ^cq) 



k=0 



(3.12) 



1 f/fj ^-l \ /fj -^-l ^^ w{z) - w{z) 

{{{J - z) eo,eo) - {{J - z) eo,eo)) = 



z — z 



z — z 



For z fixed and w = w{z) a complex parameter (|3.12| ) gives rise to an equation in the complex 
'-plane, which is in general a circle or a point. The radius of this circle is (1^^ — z\ '}2'k=Q \Pk{z) 



Proposition (3.3.6) shows that it is a point for a bounded Jacobi operator, see |l 
discussion of limit points and limit circles. 

(3.3.9) We now introduce the Green kernel for z £ 



GkM 



for a 



fl{z)Pk{z), k<l, 
fk{z)pi{z), k> I. 



Hence {Gk,iiz)}'^^Q,{Gk,iiz)}^Q G ^^(Z>o) by Proposition |(3.3.6)| , and the map 

oo 

f{Z>o) B V ^ {Gz)v, {G{z)v)^ = Y, viGk,i{z) = {v,G^^) 



1=0 



is a well-defined map. A priori it is not clear that G{z) is a bounded map, but it is densely defined, 
e.g. on D(Z>o). The next proposition states that G is bounded. 



Proposition. The resolvent of J is given by {J — z) "^ = G{z) for z G 
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Proof. Since z G C\M C p{J) we know that {J — z) ""^ is a bounded operator, so it suffices to check 
(J — z)G{z) = 1^2(2^|j-) on a dense subspace 2?(Z>o) of ^^(Z>o). Now 



{{J - z)G{z)v) j^ = ^vi{akGk+i,i{z) + (6^ - z)Gk,i{z) + ak-iGk-i,i{z)) 

1=0 
k~l 

= '^vipi{z){akfk+iiz) + (bk - z)fk{z) + afc-i/fc-i(2;)) 



1=0 



+ X] f^ifi{z){'ikPk+i{z) + {hk - z)pk{z) + ak-iPk^i{z)) 

l=k+l 

+ Vk{akfk+i{z)pk{z) + {hk - z)fk{z)pk{z) + afc_i/fc(2)pfc-i(2;)) 
= Vk{akfk+i{z)Pk{z) + {bk - z)fk{z)pk{z) + ak^ifk{z)Pk~i{z)) 
= Vkak{fk+i{z)Pk{z) - fk{z)Pk+i{z)) 

using that fk{z) and Pk{z) are solutions to J f = zf for k > 1. The term at the right hand side 
involves the Wronskian, see |(3.1.5)| . Since \p,p] = 0, and [p,r] = —1 and fkiz) = w{z)pk{z) + rk{z), 
it follows that [f,p] = ak{fk+i{z)pk{z) - fk{z)Pk+i{z)) =1. □ 

(3.3.10) We can now find the spectral measure of J from the resolvent, see the Stieltjes-Perron 
inversion formula (2.2.4); 

-E-Mtjffa, 6)) = limlim / (Gix + ie)u,v) — (Gix — ie)u,v) dx. 

' ^ ^ <5iO e[0 2m J^^s 

First observe that 

oo 

{G{z)u,v) = ^ Gk,i{z)uiVk = ^ fi{z)pk{z)uiVk + ^ fk{z)pi{z)uiVk 

k<l k>l 

= ^ fi{z)pk{z){uiVk + UkVi){l - -6k,i), 



k,l=0 



k<l 



by splitting the sum and renaming summation variables. The factor {l—2^k,l) is introduced in order 
to avoid doubling in the case k = I. Since fi{z) = w{z)pi{z) +ri{z)., withp; and n polynomials, the 
only term contributing to {G{x + ie)u,v) — {G{x — ie)u,v) as e | comes from w^z), cf. (3.1.2) 
Hence, 

Eu,v{{a,b)) = 

limlim-— / [w{x^ie)-w{x-iE)\^pi(x)pk{x){uiVk^UkV\)(\--bk,i)dx 
m eio 2-Ki J^^s ^ ^ ^ 2 

and using Proposition (3.1.2) and symmetrising the sum gives 

Eu4{a,b))= [ {Uu}{x){Uv){x)d^x{x), 



(a,b) 



(3.13) 



where U: £^(Z>o) -^ L'^il^) is the unitary operator U of Theorem (3.3.5) . Note that ( p. 13 ) proves 
that f/ is a unitary operator by letting a — > — oo and 6 — > oo, so that 



{u,v) = / {Uu){x){Uv){x)dfi{x) = {Uu,Uv)l2(^^). 
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We can think of U as the Fourier transform; 

Uk= I {Uu){x)pkix)dfi{x), 

JR 

SO that u £ £^(Z>o) is expanded in terms of (generahsed) eigenvectors {pfc(a^)}fcLo °f ^^^ Jacobi 

QDGrcitor J 

Theorem. (3.3.11) There is one-to-one correspondence between bounded Jacobi operators and 

probability measures on M with compact support. 

Proof. Since a probabihty measure fi of compact support has finite moments, we can build the 
corresponding orthonormal polynomials and write down the three-term recurrence relation by The- 
orem (3.1.3). This gives a map /U i— > J. Using the moment generating function, analytic at cxd, and 



the Stieltjes inversion formula, cf. proof of Theorem |(3.3.5)i we see that the map is injective. In 
Theorem (3.3.5) the inverse map has been constructed. D 

3.4. Unbounded Jacobi operators. 

(3.4.1) We now no longer assume that J is bounded. This occurs for example in the following 
case. 

Lemma. Let n be a probablity measure with finite moments. Consider the three-term recurrence 
relation of Theorem (3.1. 3)| and the corresponding densely defined Jacobi operator J. If supp(/i) is 
unbounded, then J is unbounded. 

Proof. Suppose not, so we assume that J is bounded. Using the Cauchy-Schwarz inequality and 
(3.11) we obtain 

||jf'>KJ''eo,eo)| = | [ x^^d^i{x)\ 

„2fc 



This implies 



> / x'^dfiix) > (|J|| + 1)^'' f dfi{x). 

|x|>||J||+l •^|a;|>||J||+l 



lUII \'^'' f 
\J\\ + ^'' Jla'l>IIJII+l 



and by the assumption on fi the right hand side is strictly positive. The left hand side tends to 
as /c — > oo, so that we obtain the required contradiction. D 

(3.4.2) We use the notions as recalled in §2.3. Recah from ( |3.10| ) that J is a densely defined 
symmetric operator. Let us extend J to an arbitrary vector v = YlT=Q'^k&k G ^^(^>o) by formally 
putting 

oo 

J* v = {aovi + bovo) eo + ^(afc^fc+i + hvk + au-iVk-i) e^, (3.14) 

fc=i 

which, in general, is not an element of £^(Z>o). 

Proposition. The adjoint of {J,V{'L>q)) is {J*,V*), where 

V* = {ve ^^(Z>o) I J*v G f{Z>o)} 

and J* is given by (3.14). 



Proposition (3.4.2) says that the adjoint of J is its natural extension to its maximal domain. In 
case J is essentially self-adjoint we have that (J*, 2?*) is self-adjoint and that (J*,P*) is the closure 
of J. 
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Proof. To determine the domain of J* we have to consider for which v G £^(Z>o) the map w 
{Jw,v), D(Z>o) ^ C is continuous, see |(2.3.T) . Now, with the convention e_i = 0, 

\{Jw,v)\ = I ^u^fc(«fcefc+i + bkek + ak-iek-i,v)\ 

k 

= I '^Wk{akVk+i + bkVk + ak-iVk-i)\ 



k 



< \\w\\w^MkVk+i + bkVk +aA:--i^'fc-in^ = ||it'||||J*w| 



k 



1 



where the sums are finite since w £ T>{'L-:>q) and where we have used the Cauchy-Schwarz inequahty. 
This proves w i-^ {Jw,v) is continuous for v G D*. Hence, we have proved D* C T>{J*). 
For the other inclusion, we observe that for v G T){J*) we have 



y^ Wk{akVk+i + hvk + ak-iVk-i) \ < C\\ 



_ }w\ 

k 

for some constant C independent of u; G P(Z>o). Speciahse to Wk = iJv)k iov < k < N and 

Wk = ioi k > N, to find \\w\\'^ < C||tt;||, or 

N 



(E 



a^Vk+i + hkVk + ak-iVk-i\^Y < C 



k=0 

and since C is independent of w, C is independent of N. Letting A^ ^ oo we find 11^*^11 < C, or 
vGV*. D 

Since the coefficients are real, J* commutes with complex conjugation. This implies that the 
deficiency indices are equal. Since a solution of J* f = z f is completely determined by /o = (/, cq), 
it follows that the deficiency spaces are either zero-dimensional or one-dimensional. In the first case 
we see that J with domain P(Z>o) is essentially self-adjoint, and in the second case there exists a 
one-parameter family of self-adjoint extensions of (J, P(Z>o)). Since the only possible element in 
Nz is {Pk{z)}'^Q we obtain the following corollary to Proposition |(2.3i4) . 



Corollary. (J, I?(Z>o)) is essentially self-adjoint if and only ifYl'k^=o IPfc(-^)P ~ "-^ f^''" ^^^ ■^ ^ 
if and only if YlT=o \Pk{z)\'^ = oo for some z G C\IR. 

(3.4.3) In case J has deficiency indices (1,1) every self-adjoint extension, say (Ji,2?i), satisfies 
(J,P(Z>o)) C (Ji,Pi) C (J*,P*), and by Proposition |(3.4.2)| (Ji,Pi) is the restriction of J* to a 
domain Pi. Since the polynomials pk and r^ are generated from the restriction to DCljyo), these 
polynomials are independent of the choice of the self-adjoint extension. To describe the domains of 
the self-adjoint extensions we rewrite the sesquilinear form B of (2.3.4) in terms of the Wronskian; 
using the convention n_i = = v_i 
N 



^{J*u)kVk - Uk{J*v)k 
k=0 

N 



^{a-kUk+i + bkUk + ak-iUk-i)vk - Uk{akVk+i + bkVk + ak-iVk-i) 



k=0 



N 



= [u,v\o + 2_^\^[u,v\k - [U,V\k-^lj = [u,v\n, 

k=l 

where [n, uj^ = ak{uk+iVk — Vk+iUk) is the Wronskian, cf. (|3.7]). So for u,v G V* we have 
B{u,v) = liiaiN^oc[u,v]N- Note that the limit exists, since u,v & V* . 
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Lemma. Assume that J has deficiency indices (1,1), then the self-adjoint extensions are in one- 
to-one correspondence with {J*,'Dg), 6 £ [0, 27r), where 

Ve = {v£V*\ lim [v, e^Vi + e-^>-i]7v = 0} 

A''-^oo 



where J*'ip±i = ±iil^±i, ii'i)k = {ip-i)k and HV'iill = 1- 

So (V±i)fc = Cpk{±i) with C7-2 = J2Zo \Pk(i)\^ < oo by Corollary pX2]| . 

Proof. Since n+ = n_ = 1, Proposition |(2.3.4)| gives that the domain of a self-adjoint extension is 
of the form u + c(V'i + e-2*V-i) with u£ V{J**), c G C. Using B(V'i,V'i) = 2i, B{ip^i,ip_i) = -2i 
and B{ipi,ip-i) = shows that 

B{u + c(Vi + e-2^V_,), e^Vi + e"'V-i) = 5(^^, e^Vi + e-'U-i) 
= {ru, e^V. + e-*V_i) - (n, e^V>i + e"^^ J>_i) . 

Now observe that for v G A'^-i-j we have B{u, v) = zizi{u, v)j* = 0, u £ P( J**), using the graph norm 
as in Proposition |(2.3.4)| . D 



The condition {ipi)k = {ip-i)k in Lemma K3.4.3)| is meaningful, since J* commutes with complex 
conjuagtion. It is needed to ensure the one-to-one correspondence. 

(3.4.4) The results of §3.3 go through in the setting of unbounded operators up to minor changes. 
Lemma K3.3.2) remains valid, as it has nothing to do with the unboundedness of J. Lemma |(3.3.3)| 
implies that for an unbounded Jacobi matrix we have to deal with unbounded sequences {a^}, 
{6fc}, or at least one of them is unbounded. (3.3.4) goes through, using Theorem (2.4. 1)| , after 
choosing a self-adjoint extension of J. This operator is uniquely determined if {n+,n-) = (0,0) 
and is labeled by one real parameter if (n-|-,n_) = (1,1). For this self-adjoint extension |(3.3.4)| 
and Theorem (3.3.5) remain valid, except that ^ is no longer compactly supported. In this case we 



observe that /i has finite moments, since jjg x dfi{x) = {J cq, cq) < oo. From the unicity statement 



in Theorem (2.4.1) we deduce that the polynomials are dense in L^(^), which is needed to obtain 
the unitarity of U in Theorem (3.3.5) . The results of §3.3 until Theorem [(3.3.11) remain valid after 
replacing J by a self-adjoint extension. Observe that / as in Proposition (3.3.6) is contained in 



the domain of a self-adjoint extension, so that the calculation in Proposition (3.3.9) remains valid. 



In the unbounded case it is no longer true that ( 3.12| ) in the complex t^-plane describes a point. 
Theorem |(3.4.5)| and Corollary (3.4.2) show that in the general case ( 3.12D describes a circle (with 
positive radius) if and only if the Jacobi operator has deficiency indices (1,1) if and only if the 



corresponding (Hamburger) moment problem is indeterminate. Finally, Theorem (3.3.11) cannot 
be extended to unbounded self-adjoint Jacobi operators. 

(3.4.5) There is a nice direct link between the moment problem in §3.2 and the Jacobi operators. 
It gives an answer to question 2 of (3.2.1). 



Theorem. The (Hamburger) moment problem is determinate if and only if the corresponding 
Jacobi operator is essentially self-adjoint. 

Proof. Assume first that the corresponding Jacobi operator is essentially self-adjoint. Then (J — 
z)(P(Z>o)) is dense in i^{Z>o) by Proposition |(2.3.4)| , z £ C\M. The density and Lemma [(3^3^ 
say that we can find polynomials Rk{x] z) such that ||(J — z)Rk{J] z)eo — eo|| — > as A; — > oo. Let 
H be any solution to the moment problem, so that J is realised as the multiplication operator on 
the closure of the span of the polynomials in L'^[fi). Then 



lim 



\{x — z)Rk{x; z) — 1| dfi{x) = 0. 
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Since z G C\M we see that x 



-1 ■ 



is bounded for x G M, hence 



fc— >oo 



\Rk(x\z) 



1 



X — z 



' dfi{x) = 0. 



Since /i is a probabihty measure L'^i^J-) ^^ L^il^), so the Stieltjes transform of fi is given by 

= hm / Rk[x\ z) dii{x) , 

X-Z k~>ooj^ 

but the right hand side only involves integration of polynomials, hence only the moments, so it 
is independent of the choice of the solution of the moment problem. By the Stieltjes inversion of 
Proposition (3.1.2) this determines /j. uniquely. 

For the converse statement, we assume that J has deficiency indices (1, 1). Let Ji and J2 be two 
different self-adjoint extensions of (J, P(Z>o)). We have to show that they give rise to different 



solutions of the moment problem. By the Stieltjes inversion of Proposition (3.1.2) it suffices to 



show that the Stieltjes transforms ((Ji — z)~^eo,eo) and ((J2 — z)~^eo,eo) are different. 

Let us first observe that cq ^ Ran(J**— z). Indeed, suppose not, then there exists u £ V{J**) such 
that (J** — z)u = eo, and taking the inner product with d ^ v £ N^ gives {eo,v) = {{J** — z)u, v) = 
{u, (J* — z)v) = 0, since u £ 'D{J**) and v £ 2?(J*) and v £ N^. So v is an eigenvector of J* with 
1^0 = 0, so that V = by (3.4.2)| . This contradicts the fact that J has deficiency indices (1,1). 

1,2, by 



From the previous paragraph we conclude that (Jj — z) cq £ 'D{J*)\D{J** 



Proposition |(2.3.4)| . Now the dmi[V{Ji)\V{J**)) 
T>{Ji) = V{J2) and hence Ji 



1, so that (Ji 



eo 



{Ji — z) -"^eo implies 
J2 by Proposition |(2.3.4)| . So the vectors ( Jj — z)~^eo are different. 
To finish the proof we need to show that the zero-components are also different. Let u = {Ji — 
^eo - (J-?. - zy 



z)-^' 



(J2 — z) ^eo / 0, then 

{J* - z)u = {J* - z){Ji 



) eo - {J* - z){J2 - z) eo = eo-eo = 0, 



or u e Nz- By (3.4.2)| {u,eo) = implies u = 0, so (n, eo) / 0. Hence, Ji and J2 give rise to two 
different Stieltjes transforms, hence to two different solutions of the moment problem. D 



(3.4.6) From Lemma (3.3.3) we see that boundedness of the coefficients in the three-term re- 
currence relation of Theorem (3.1.3) for the orthonormal polynomials implies that J extends to 
a bounded self-adjoint operator. For unbounded coefficients there are several conditions on the 
sequences {a^}, {6^} that ensure that the Jacobi operator J is essentially self-adjoint, and hence, 
by Theorem (3.4.5), the corresponding moment problem is determinate. We give three examples. 



Proposition, (i) // X^^q IT ~ ^'°' ^^e'^ J ^^ essentially self-adjoint. 

(ii) // aj. + hk + Ofc-i < M < CO for k £ 'N, or if a^ — h^ + Ofc-i < M < co for k £'H, then J is 

essentially self-adjoint. 

(iii) If Yl'k^=o''^2k ~ °*^' '"^^fi'^e TTij is the j-th moment trij = J^x^ dfi{x), then J is essentially 
self-adjoint. 

Proof. To prove (i), we use the Christoffel-Darboux formula (l3.8| ) with x = z, y = z, and using 
that the polynomials have real coefficients, we find for z £ 



N 



k=0 



k[Z) 



ON 



{pn+i{z)pn{z) -pn{z)pn+i{z)). 



Since the left hand side is at least 1 = po{z), we obtain 



1 < 



UN 



■{\pn+i{z)\\pn{z)\ + \pn{z)\\pn+i{z)\) 



2aN 



\z — z\ 



'\PNiz)\\pN+liz)\. 
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For non-real z we have C = 2l\z — z\ > and so 

oo oo oo 

Y.- <CY^\v^{z)\W^^{z)\<cY^\v^{z)\^ (3.15) 

fc=0 ^ k=0 k=0 

by the Cauchy-Schwarz inequahty (2.1.1) in £^(Z>o). So, if YlT=o i~ diverges, we see from Corollary 



(3.4.2) and ( 3.15| ) that J is essentially self-adjoint. 



The cases in (ii) are equivalent, by switching to the orthonornial polynomials (— l)"'p„(— x) for 
the measure fj.{A) = fi{—A) for any Borel set A C M. We start with an expression for Pn{x) in 
terms of the lower degree orthonornial polynomials; 



ra-l 



Pn{x) = 1 + X] — Y,i^ ~ ^j ~ «i ~ "■j-i)Pjix)^ (3-16) 

with the convention a„i = and p^i(x) = 0. In order to prove ( |3.16| ) we start with the telescoping 



series 



n-l 



an-i{pn{x) -pn-i{x)) = ^ (aj(pj+i(x) - Pj{x)) - aj-i{pj{x) -pj-i{x))) 



j=0 
n-l 



■ y^{x — bj — aj — fflj-i) Pj{x) 

j=0 



by (0), (|33). This gives 

Pn{x) = Pn-l{x) H ^(X - hj - aj - aj-i)pj{2 



n-l 



O-n-l 



j=0 



and iterating this expression gives ( 3.1£| ). 

Now take x > M in (3.16) to obtain inductively Po(^) = 1) Pk{x) > 1 for A: G N, so that 
SfcLolPfc(^)P i^ divergent. Since the polynomials pk have real zeroes, see Lemma (3.3.7) , we 
have pk{x) = CkY\i^i{x — xi) with Xj G M which implies that \pk{x + iy)\ > \pk{x)\. Hence, 
YlT=o \Pk{x + iy)]"^ = oo and Corollary (3.4.2) shows J is essentially self-adjoint. 



Finally for (iii) we note that 

1 = WPkWmp) < \\HPk)x''\\ = HPk)yM2k 



^/m2k 



aoaia2 • . . ak-i 
by the triangle inequality, where \c{pk) denotes the leading coefficient of the polynomial pk- Hence 

oo k—l 



m. 



-l/2k 
2k 



<rlll 



1 sl/fc 



ao ai 02 flfc-i' 
Using the geometric-arithmetic mean inequality, 



E 

fc=0 



-l/2k . 
™2fc < 



E(n^)" 

fc=0 J=0 



'Wa\' <Y^Piai, ^pj = l, 



1=1 



i=\ 



t=\ 



we find 



n^)""-('o-'"(n^)"^(.o-'"^E 



4 = 



1 ^=-1 ■ I 1 

1 V^ « + 1 



j=0 



i=0 
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Next use k^ < e^kl, which can be proved by induction on k and (1 + ^)^ < e, to see that we can 
estimate 

oo oo fc— 1 . ^ oo ^ oo . ^ 

fc=0 fc=0 i=0 i=0 fc=i 



and the inner sum over k can be estimated independent of i. Now the resuh fohows from (i). D 

(3.4.7) A famous example of an indeterminate moment problem has been given by Stieltjes in his 
posthumously published memoir [21|. Consider for 7 > — Stieltjes considered the case 7 = 1 — 

\ x"e-^^" "sin(27r72lnx)dx = e^^ / e-^^'"""^) sin(27r72 Inx) — 
JO JO ^ 

and put y = 7lnx — ^y— , x^^dx = "f^^dy to see that the integral equals 

J ^e 47^ / e ^ sin(27r7y + 7r(n + l))(i2/ = 0, 
J —00 

since the integrand is odd. In the same way we can calculate the same integral without the sine-term 
and we find that for any — 1 < r < 1 we have 

/" x" e-^' 1"' ^ (1 + r sin(27r72 In x)) dx = ^ exp (^^^i^) , 

so that the probablity measure 

dfi{x) = -^e^ e-^' ''^' ^ dx, X > 0, 



corresponds to an indeterminate moment problem. The corresponding orthogonal polynomials 
are known as the Stieltjes- Wigert polynomials. Put q = exp(— 1/27^) or 7"-^ = —2 In (7, then 
the orthonormal Stieltjes- Wigert polynomials satisfy ( |3.4| ), ( |3.5D with a^ = q^^ ^^ \J\ — g'^+^ and 
6fc = g~^'^(l+g~^ — g'^), see e.g. [^, §VI.2]. Note that < g < 1, and that a^ and 6^ are exponentially 
increasing. 

(3.4.8) Another example is to consider the following measure on R 

d^(x) = Ca^-y e"''''^' dx, a, 7>0. 

The moments and the explicit value for Co^-y can be calculated using the F-function; 

00 

x^-^e-''^ dx = 6"T(c), c> 0, K6 > 0. (3.17) 



For < a < 1 the (Hamburger) moment problem is indeterminate as we can see from the following 
observations. Note that all odd moments vanish. Put x = y" and c = (2n + l)/a to find from 
( 3.17] ), after doubling the interval, 

2n -felul" r 2n+l . 2n + 1 X 



By taking the real parts we obtain 

« r y2„^-«%r,os(-96|yr)dy = F(^^^±^)|6r^cos(-arg(6)^^^±^), 

so that the right hand side is zero for all n G Z>o if arg(6) = ^air. Since we need ^b > this is 
possible if |a| < 1. So for < a < 1 the (Hamburger) moment problem is indeterminate, since we 
obtain more solutions in the same way as for the Stieltjes- Wigert polynomials in (3.4.7). 
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The (Hamburger) moment problem is determinate for a > 1, which fohows using Proposition 



(3.4.6) (iii). See also Deift JS], p. 34] for another proof, in which the essentially self-adjointness of 
the corresponding Jacobi operator is established. 

Restricting the measure d^{x) to [0, oo) gives a Stieltjes moment problem, see §3.2, and for the 
Stieltjes moment problem this is an indeterminate moment problem if and only if < a < 21 see 
||19| , Ch. 1, §8]. In particular, for 2 — ck < 1; the Hamburger moment problem is indeterminate, 
whereas the Stieltjes moment problem is determinate. In terms of Jacobi operators, this means 
that the corresponding Jacobi operator J has deficiency indices (1,1), but that J has a unique 
positive self-adjoint extension, see Simon [pO| . 

(3.4.9) A great number of results exist for the description of the solutions to an indeterminate 
moment problem. A nice way to describe all solutions is due to Nevannlina; see e.g. [Q], ||2^. There 
exist four entire functions A, B, C and D, that can be described using the polynomials pk and r^, 
such that the Stieltjes transform of any solution n to the indeterminate moment problem is given 

by 

d/i(x) A{z)c^{z) - Ciz) 



x-z B{z)(j){z) - D{z) 

where <f) is any function that is analytic in the upper half plane having non-negative imaginary 
part, or = 00. 

Now the solutions to the moment problem corresponding to <J3{z) = i G MU{oo} are precisely the 
measures that can be obtained from the spectral measure of a self-adjoint extension of the Jacobi 
operator J. In particular, this implies that these orthogonality measures are discrete with support 
at the zeroes of an entire function B{z)t — D{z), so that there is no accumulation point. According 
to a theorem of M. Riesz (1923) these are precisely the measures for which the polynomials are 
dense in the corresponding L^-space, cf. Theorem (3.3.5) for the unbounded case. See Akhiezer |l|] 
and Simon [^1 for more information. 
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4. Doubly infinite Jacobi operators 
4.1. Doubly infinite Jacobi operators. 



(4.1.1) We consider an operator on the Hilbert space ^ (Z), see (2.1.2) , of the form 



Lek = ak e/c+i + bk e^ + at-i ek-i, ak > 0, 6^ € 
where {ek}kez is the standard orthonormal basis of £'^(Z) as in (2.1.2). If a, = for some i £ 



then L splits as the direct sum of two Jacobi operators as defined in (3.3.1). We call L a Jacobi 
operator on £^(Z) or a doubly infinite Jacobi operator. 

The domain ^^{L) of L is the dense subspace 2?(Z) of finite linear combinations of the basis 
elements e^. This makes L a densely defined symmetric operator. 
(4.1.2) We extend the action of L to an arbitrary vector v = Ylh=-oo '^k^k £ ^^(Z) by 

oo 

L*v= y^ {akVk+i + bkVk + ak^iVk-i)ek, 

k=—oo 

which is not an element of £'^{Z) in general. Define 

V* = {ve f{Z) I L*v G e^{z)}. 

Lemma. {L*,V*) is the adjoint of{L,V{Z)). 



The proof of this Lemma is the same as the proof of Proposition (3.4.2) . 
(4.1.3) In particular, L* commutes with complex conjugation, so its deficiency indices are equal. 
The solution space of L*v = zv is two-dimensional, since v is completely determined by any initial 
data {vn-i,Vn) for any fixed n £ Z. So the deficiency indices are equal to {i,i) with i S {0,1,2}. 
From §2.3 we conclude that L has self-adjoint extensions. 

4.2. Relation to Jacobi operators. 

(4.2.1) To the operator L we associate two Jacobi operators J^ and J~ acting on £^(Z>o) with 
orthonormal basis denoted by {/fe}fceZ>o i^ order to avoid confusion. Define 



J^fk 



J^fk 



akfk+i + hfk + ak-ifk-i, for A; > 1, 
^ao/i + 6o/o, for A; = 0, 

a-k-2 fk+i + b-k-i fk + a-k-i fk-i, for k > 1, 
a~2 /i + ^-1 /o, for k = 0, 



and extend by linearity to P(Z>o), the space of finite linear combinations of the basis vectors 
{/fclfcLo of £^(Z>o). Then J^ are densely defined symmetric operators with deficiency indices (0,0) 
or (1, 1) corresponding to whether the associated Hamburger moment problems is determinate or 



indeterminate, see §3. The following theorem, due to Masson and Repka [17|, relates the deficiency 

Tneorem. (4r.2."2) '(Masson and Repka) The deficiency indices of L are obtained by summing the 
deficiency indices of J^ and the deficiency indices of J~ . 

Proof. Let Pk{z), Qk{z) be generated by the recursion 

zpkiz) = akPk+i{z) +bkPk{z) +ak-iPk-i{z), (4.1) 

subject to the initial conditions Pq{z) = 1, P-i{z) = and Qq{z) = 0, Q-i = 1. Then, see 



(4.1.3) , any solution of L*v = zv can be written as Vk = VQPk{z) + v-.iQk{z). Note that {p^}°° 



A:=0' 



p].{z) = Pk{z), are the orthonormal polynomials corresponding to the Jacobi operator J^ and that 
{Pfc IfcLo' Pk(-^) ~ Q-k-i{z), are the orthonormal polynomials corresponding to J~ . 
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Introduce the spaces 



-1 



A:=— oo 

oo 

St = {{Mr=-oo \L*f = zf and J] \h\' < oo}. 



(4.2) 



fc=0 



From (4.1.3) we find dimSj < 2. The deficiency space for L is precisely Sf n S~ . 

From the results in §3, in particular Corollary (3.1.4) , Proposition (3.3.6) , Corollary (3.4.2)| and 
(3.4.4) , we have for z G C\M that dimS*^ = 1 if and only if J^ has deficiency indices (0,0) and 
dimS'j^ = 2 if and only if J has deficiency indices (1, 1). 

Let us now consider the four possible cases. Case 1: J"^ and J" have deficiency indices (1,1). 
Then dim 5^ = 2 and by K4.1.3)| we get S^ = S~ , so that L has deficiency indices (2,2). Case 
2: one of J is essentially self-adjoint. We can assume that J'^ has deficiency indices (0, 0) and 
has deficiency indices (1, 1). Then dim 5*+ = 1, dim5~ = 2 and by (2.1.6) S~ coincides 

zv, i.e. every solution of L*v = zv is square summable at — oo. 
So (4^) shows that Sf C S~ and hence dimS*^ n 5~ = 1. So L has deficiency indices (1, 1) in 
this case. Case 3: J^ and J^ have deficiency indices (0,0). Then diuiSt = 1 and we have to 
show that 5"+ Ci S^ = {0}. Let v^ = {v^, }'^^_^ span the space Sf, and we have to show that 

v~^ is not a multiple of v~ . Let Vj, = C^Pk{z) + C^Qk{z). We calculate C^ and C^ explicitly 
in terms of the Stieltjes transform w^{z) = f^{x — z)~^ dfj,^{x) of the spectral measures /i^ for 
J^. Since J~^ and J~ are both essentially self-adjoint, the spaces Sf are described in Proposition 



that J 

with the solution space of L*v 



(3.3.6) . It follows from the recursion (4T) that the associated polynomials r^ (z) for J^, see (3.1.4) , 



satisfy r'^iz) 

4 = f^i^) - 

w~{z)Qk{z)- 

91 = 



-1 



Qkiz), k > 0, and r^ 



-1; 



-a_{P_ 



k-i[z) 



k > 0. By Proposition |(3.3.6)| 

w+{z)Pk{z) - a^^Qkiz), k>0, and w^ = flk-i(z) = W iz)pZk-i{z) + rZk-ii^) = 
aZiPk{z). So d^l = w^{z), Cjf. = ^, C^ = -^ and C_ = w^{z), and consequently 



-Ol w 



0(1), 



oo, 



C^ 

C^ 



a^iw {z) 



0[z\ 



oo. 



In case v^ is a non-zero multiple of v the quotients have to be equal. Since they are also inde- 
pendent of z we find the required contradiction. This implies S^ H S^ = {0} or L has deficiency 
indices (0,0). D 



(4.2.3) As in [(3.4.3) we can describe the sesquilinear form B as introduced in (2.3.4) in terms of 
the Wronskian [u, f]^ = o-kk^k+x'^k — UkVk+i)- Note that, as in §3, the Wronskian [u,v] = [u,v]k is 
independent of k for L*u = zu, L*v = zv, and then [u,v] 7^ if and only if u and v are linearly 



independent solutions. Now, as in (3.4.3 

N 



k=M 



y^ {L*u)kVk - Uk{L*v)k 

VI 

N 

y^ {akUk+i + Kuk + ak-iUk-i)vk - Uk{akVk+i + bk^k + ak-iVk-i] 



k=M 

N 

^ [u,v\k - [u,v\k-i = [u,v\n - [u,v]m-i, 

k=M 
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SO that 

B{u,v) = lim [u,v]n — lim [u, -uJm, n, uGP*. 

In particular, if J^ is essentially self-adjoint, we have liinM^-oo[u, v]m = 0, so that in this case the 
sesquilinear form B is as in §3. In case J"*" has deficiency indices (1,1) we see as in Lemma (3.4.3) 
that the same formula for the domains of self-adjoint extensions of L are valid. 

Lemma. Let J~ , respectively J^ , have deficiency indices (0,0), respectively (1,1), so that L has 
deficiency indices (1,1). Then the self-adjoint extensions of L are given by {L*,'Dg), S [0, 27r), 
with 

Vg = {veV*\ lim [v, e'^^i + e-'^^-i]N = 0} 

A''— >oo 



where L*^±i = ±i^±i, {^i)k = ($-i)fc (^^d \\^±i\\ = 1. 

The proof of Lemma |(4.2.3)| mimics the proof of Lemma K3.4.3) 

4.3. The Green kerneL 

(4.3.1) From on we assume that J~ has deficiency indices (0,0), so that J~ is essentially self- 



adjoint and by Theorem (4.2.2) the deficiency indices of L are (0,0) or (1, 1). (The reason for the 
restriction to this case is that in case L has deficiency indices (2, 2) the restriction of the domain 
of a self-adjoint extension of L to the Jacobi operator J^ does not in general correspond to a 
self-adjoint extension of J^, cf. [0, Thm. XII. 4. 31].) Let z G C\M and choose $2 G 5^, so that $^ 



is determined up to a constant. We assume {^z)k = (^z)fc, cf. Lemma (4.2.3) 



(4.3.2) Let (fiz G S'^ , such that {4>z)k = {(f>z)k- We now show that we may assume 

1. [^z,^z]^0, 

2. (pz, defined by {(j)z)k = for k < and {(j)z)k = {4'z)k for A; > 0, is contained in the domain of 
a self- adjoint extension of L. 

First observe L*4)z = z<f)z + a_i(((/)2)oe_i — ((^2)_ieo), so that 0^ G T>* . In case L is essentially 
self-adjoint, (L* ,T>*) is the unique self-adjoint extension and (2) is valid. In this case (1) follows 
from case 3 in the proof of Theorem K4.2.2)| . 

In case L has deficiency indices (1,1), we have S~ C S'^ and (2) implies (1). Indeed, if (2) 
holds and (1) not, then $2 G £^(Z) and ^z = C(f)z is in the domain of a self-adjoint extension. 
Since L*^z = z^z this shows that the self-adjoint extension would have a non-real eigenvalue; a 
contradiction. To show that we can assume (2) we observe that 

lim [4>z,e'^^i + e-'^<^^i]N= lim [</>„ e'^^i + e'^^^.i]^. 



which exists since (j)z, ^±i ^T>* . If the limit is non-zero, say K^ we use that 

-i\N 

-id 



A{z,e)= lim [$2,e^^$i + e-^^$_i 



= {L*<^z,e''^i + e-''<^-,) - {^z,L*{e''^i + e-''<^^,)) 

= e-'\z + i)($„ ^i) + e\z - i)($„ $_i) / 0. 

Otherwise, as before, $2 would be in the domain of a self-adjoint extension of L by Lemma (4.2.3) 
So that replacing <^z by (\>z — ~Mrff)^^ gives the desired result, since S~ C S^ . 
(4.3.3) Let {L' ,T)') be a self-adjoint extension of L, assuming, as before, that J~ has deficiency 



indices (0, 0). Let t/i^ G 5+, $2 G S^ as in |(4.3.2)| . We define the Green kernel for z G C\M by 

^ . .^_}^U^z)k{<Pz)i, k<l, 
'^'^'' [</>., ^.] 1 ($2)/ (</'.)fc, k>l. 
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So {Gk,i{z)}^=_^AGkA^)}r=-^ G ^{1) and e{Z) 3 v ^ G{z)v given by 

oo 

{G{z)v)k= Y. viG{z)k,i = {v,G^z)) 

l=—oo 

is well-defined. For v £ P(Z) we have G{z)v £ T>' . 

Proposition. The resolvent of [L',!)') is given by {V — z)^^ = G{z) for z G C\M. 

Proof. Note C\IR C p{L'), because L' is self-adjoint. Hence {L' — z)~^ is a bounded operator 
mapping £^(Z) onto V . So u i-^ ((L' — z)~^v)k is a continuous map, hence, by the Riesz repre- 
sentation theorem, ((L' — z)~^v)k = {v,Gk^.{z)) for some Gk^.{z) G £^(Z). So it suffices to check 



{V — z)G{z)v = V for V in the dense subspace T)['L). As in the proof of Proposition (3.3.9) we have 

fc-i 



,$^]((L'-z)G(2;)t>)^= ^ v/(afc((/.^)fc+i + (6fc-z)((/.^)fc + afc_i(</)^)fc_i)($^ 



Z=— oo 



+ XI vi{ak(^^)k+i + {hk- z){^^)k + ak-i{^z)k-i){ct)z)l 



l=k+l 



+ Vk{ak{^z)k{(t>z)k+i + {bk - z){^z)k{(pz)k + afc-i(^z)fc-i(0z)fc) 

=Vkak{{'^z)k{<Pz)k+l - {'^z)k+l{(pz)k) = Vk[4>z,'^z] 

and canceling the Wronskian gives the result. D 

(4.3.4) As in |(3.3.10)| we can calculate 



°° 1 1 

{G{z)u,v) = ^ Gk,i{z)uiVk = . , y^X^z)k{4>z)i{uiVk + UkVi){l - -Sk,i), (4.3) 






but in general we cannot pin down the terms of ([4.3|) th at will contribute to the spectral measure 
of V in the Stieltjes- Perron inversion formula of |(2.2.4)| , (2.4.1) . We now consider some examples. 



4.4. Example: the Meixner functions. This subsection is based on Masson and Repka |17|, 
see also Masson [|^] . We extend to the results of [|^] by calculating explicitly the spectral measure 
of the Jacobi operator. 



(4.4.1) In this example we take for the coefficients of L in (4.1.1)| the following; 



Ofc = ak{a, A, e) = \/(A + k + £ + l){k + e - A), hk = bk{a, A, e) = 2a{k + e). 

Without loss of generality we may assume that < e < 1 and a > by changing to the orthonormal 
basis fk = {—l)^ek and the operator —L. We will do not so in order to retain symmetry properties. 
The conditions a^ > and 6fc € M are met if we require a G M, e € M for 6^ E R and either 
X = -^ + ib, b > or e G [0, i) and -i < A < -e or e e (i 1) and A G (-A,e - 1) for a^ > 0. 



Now L is essentially self-adjoint by Theorem (4.2.2) , since J are essentially self-adjoint by (i) of 



Proposition (3.4.6) 



(4.4.2) In case we choose A, e in such a way that Oj = for some i G Z, the corresponding Jacobi 
operator on the C-span of Ck, k > i, can be identified with the three-term recurrence relation for 
the Meixner-Pollaczek, Meixner or Laguerre polynomials depending on the size of a. In case aj ^ 0, 
we can still consider the corresponding Jacobi operator on ^^(Z>o) by putting a_i = 0. Then the 
corresponding orthogonal polynomials are the associated Meixner-Pollaczek, Meixner or Laguerre 
polynomials. 
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(4.4.3) First note that if we find a solution ui.{z) = Uk{z; a, A, e) to 

zuk{z) = akia,X,e)uk+i{z) + bk{a, X, e) Uk{z) + ak-i{a, X, e) Uk-i{z) (4.4) 

then Vk{z) = {—l)uk{—z) satisfies 

zvk{z) = ak{-a,X,e)vk+i{z) +bk{-a, X,e)vk{z) + ak-i{-a, X,e)vk-i{z) 

and Wk{z) = u_fc(z; — a, A, — e) also satisfies (|4.4|). Indeed, introducing a new orthonormal basis 
fk = e_fc of ^^(Z) we see that L is given by 

L fk = O-fc-l fk+l + b-k fk + d-k fk-l 

and 

a-k-iiX,e,a) = ak{X,-£,-a), b^kiX,e,a) = 5fc(A, -e, -a). 

(4.4.4) In order to find explicit solutions to L*v = zv we need the hypergeometric function. 
Recall 

{a)k = a{a + l)...{a + k-l)= ^ / , 

r(a) 

where the series is absolutely convergent for \x\ < 1. The hypergeometric function has a unique 
analytic continuation to C\[l,oo). It can also be extended non-uniquely to the cut [l,oo). See e.g. 
1^, Ch. II] for all the necessary material on hypergeometric functions. 

Lemma. Let a^ > 1. The functions 



4iz) =(±2)-(y;;^)-'=^^^5+^±£±l)I(^Ay 

^ fk + e + l + X,k + £-X 1 , 
X 2i^i , , , . , ,^ r^ — 7;7T± 



A; + e + l±z/2^42-^T'2 2^^^^^ 



and 



/-k-e + l + X,-k-e-X 1 



-k-e + l± z/2V^2-^T ' 2 ^ 2V^2T7T^ 
satisfy the recursion ([4.^). 
Proof. The hypergeometric function is a solution to the hypergeometric differential equation; 

x{l — x)y" + (c — (a + 6 + l)x)y' — aby = 0. 
Now /(x; a, 6, c) = r(c) 2-P"i(Q; ^; c; x) satisfies f'{x;a, b,c) = f{x; a + 1, 6 + 1, c + 1). Hence, 

rk = fix; k + l + e + X,k + e-X,k + l+e-y), 

x(l - x)rk+i + {k + e-y-2{k + e)x)rk - {k + e + X){k + e - X- l)rk^i = 0. 



Replace x = \ — ,\ and y = z/2^a'^ — 1. Let 

rk = {-2f{a^ - 1)^ Vr(fc + e + l + X)T{k + £- X)pk 
then pk satisfies (O) . 
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This proves the lemma for u^{z). The case 'U^(^) follows by replacing a by —a, and applying 



(4.4.3) . Replacing k by —k, a by —a and e by — e and using (4.4.3) gives the other sets of solutions 
to (0). D 

From now on we assume that a^ > 1 in order not to complicate matters. The case a^ = 1 
corresponds to the Laguerre functions, and the case < a^ < 1 corresponds to the Meixner- 



Pollaczek functions, see [^] and (4.4.11' 



(4.4.5) In order to find the asymptotic behaviour of the solutions in Lemma (4.4.4) we first use 



2F,('";^x 



(l-x) 



c~a— 6 



2F1 



c — a,c — b 



so that 



'k + e + l + X,k + e- X 1 
2-^1 I . . ,_ r^ -; - ± 



k + e + l±z/2V^^^T'2 2V^^^T 



-k-e±z/2^/IiF^| 



2yfc 



^ + 0{t)), 



oo. 



Use 



to find 



^^ = --'(1 + 0(1)), 

Tiz + b) ^ z ' 



Vr(A: + A + e + l)r(A: + e-A) 



z ^ oo, I arg(z)| < vr. 



|r(A: + l + e + 2/)| 
so that we find the asymptotic behaviour 

a 



fc-%-|(l + 0(i)), 



oo. 



2 2V^2-^T ^ ^' 



i^ 



as A: — > oo. This implies that u^(z) £ Sf if a < —1 and u [z) G Sf if a > 1. 
Similarly we obtain the asymptotic behaviour 



iU 



^ ± ^e±R./2v/;?^^-|TJR./2v/^?^^, ^ 7^^^)-Vl + 0(-)), 



as A; — > oo, so that v~ {z) G Sj if a < — 1 and f+(z) G Sj if a > 1. 

Note that we have completely determined 5^, hence 0^ and <1>2, since these spaces are one- 
dimensional. 

(4.4.6) If we reparametrise the parameter a = \(^s ^ s~'^^ ^ then 



{s, s"^} = {a + \/a2 - 1, a - \/o2 - 1}. 

Note that a?' — \= 4(5— s"^)^. In this case we can let u^(z) and t;^(z) correspond to hyper geometric 
series with only s^^-dependence. The results can be written somewhat nicer after transforming 
Lemma (4.4.4)| by 



2-f^i I \x 



(1 - x)-\F^ 



a,c — b X 



We leave this to the reader. 

(4.4.7) Note that v^{z) are linearly independent solutions to ([4.4|), because they display different 
asymptotic behaviour. Since the solution space to ( [4.4| ) is two-dimensional, we see that there exist 
constants such that 



u. 



\z) = A+{z)vt{z)+B+{z)v-^{z), u^{z) = A-{z)v+{z)+B-{z)v,{z). 
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Or stated differently, relations between hypergeometric series of argument x and 1 — x. Relations 
connecting hypergeometric series are very classical. We use 

\a + o + l — c J \ (^ J \ 2 — c 

^ _ r(a + b + i-c)r(i-c) r(a + 6 + i-c)r(c-i) 



r(a + i-c)r(6 + i-c)' r(a)r(6) 

H^— fe — e + l + A, 6i-^— fc — e — A, ci-^— A; — e + : 
2-Fi is as in ^^(2), the second one as in v~l{z)^ and the third one as in % (z). This gives, using 
y = z/2^/a? - 1, 



with a H^ — /c — e + l + A, 61— > —A; — e — A, c 1— > —k — e + l + y, xi— > ^ ."^ Then the first 

^ 2va — 1 



S-(z)=(-2)-'=(Va2-l) 



-fc r(-fc + i-£-7/) 

Vr(-A; + A - e + l)r(-A: - e - A) 



from factor of v^ {z) 



T{-k- e + l + \)T{-k- e- X) ,1 a .y-k~e A a -s-y-k-e 

' T{l-k-e-y)T{-k-e + y) ^2 " 27P^^ ^2 2^^^' 

from B 






{k + l + £-y) 

from factor of Uj^{z). 

Now B~{z) has to be independent of /c. After canceling common factors we use r(z)r(l — z) = 
vr/sinvrz, z Z, to obtain 



5-(.) = (^^""7(4(1 - a2)r ^^"^"(^ - ")) 



\/a2 — 1 4- a^ Y^sin(7r(e — A)) sin(7r(— e — A)) 

Similarly, also using T{z + 1) = zT{z), we find 

TT /Va^-l-ax!/ (4(1 -a^))^ 



^-fz) 



r(l + A-y)r(-A-y)V^/^^3T + a; ^sin(^(e - A)) sin(7r(-e - A)) ' 

(4.4.8) Next we calculate the Wronskians [v~ (z) , v^ (z)] and [u~(z),f"'"(z)]. To calculate the first 
Wronskian observe that we can take the limit /c — > 00 in 

[v~iz),v^iz)] = a-k{vZk+iiz)vlk{z) - vZkiz)vtk+iiz)) 
using the asymptotic behaviour of |(4.4.5)| . This gives, using y = z/2\/a'^ — 1, 

\v'(z),v+(z)] =a_fc(- — -^ — Y~y(- + , " Y+y 

X ({k - l)y-^a - Va'^-l)^-''k-^-y{a + Va"^ - 1)^^ 

- F-5(a - Va2 - 1)-'^(A; - l)-5-?/(a + x/a2 _ i)i-fe^ 

and taking out common factors and using that a_jfc = k{l + 0{^)) gives 



\/a2 — 1 + a\y 



■(.).„+(.)! = -2sA?3T(4(i-„')r(^° - +°)'. 

^ Va — 1 — a^ 
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From 1(4.4.7)1 it follows that, y = zl2\/a^ - 1, 



\u 



\z),v^{z)\ = B-{z)[v-{z),v+{z)] = -2^^ 



1- 



sin(7r(y-e)) 



Y^sm(7r(e — A)) sin(7r(— e — A)) 



We see that the Wronskian [u (z),f+(z)], as a function of z, has no poles and it has zeroes at 

2; = 2(e + /)Va2-l, I G Z. 

(4.4.9) With all these preparations we can calculate the spectral measure for the doubly infinite 
Jacobi operator L with coefficients as in (4.4.1), where we assume a > 1. So we can take (j)^ = u~{z) 
and <1>2 = v'^{z). And we see that these solutions are analytic in z, since V{c)~^ 2Fi{a,h;c;z) is 
analytic in c. Hence, the only contribution in the spectral measure, cf. (4.3.4)| , comes from the 
zeroes of the Wronskian, so that the spectrum is purely discrete. 

Theorem. The operator L on ^^(Z) defined by (4.1.1) , (4.4. 1)| is essentially self-adjoint. For a > 1 
its self-adjoint extension has completely discrete spectrum {2(e + l)\/a'^ — l}i^i, and in particular 
the set 

{u-{2{e + 0\/a2-l) I / G Z} 

constitutes a complete orthogonal basis of i'^{'L). Moreover, 



(2(e + /)Va2-l)||2 



4(a2 - l)(a - V^?^^)Y(a - V? 



+ ^/^F^l 



1\' 



, + y^^^T 



r(e + /-A)r(l + e + / + A). 



Remark (i) Note that hyper geometric expression for u^ as in Lemma (4.4.4) displays Bessel co- 
efficient behaviour, and we can think of the orthogonality relations as a natural extension of the 
Hansen-Lommel orthogonality relations Yl'kL-oo Jk+n{z)Jk{z) = (5n,o for the Bessel functions, see 



(4.4.10) 



(ii) Note that the spectrum is independent of A. 

Proof. Put <f)z = u~{z) and ^^ = v~^{z), then we see that, see ( [4. 3D , 

°° 1 1 

^ Gk,i{z)uiVk = r , ^ . ^{^z)k{(t>z)i{uivk + UkVi){l - -6k,i) 



{G{z)u,v) = 

k,l = — OD 

Hence, the only contribution in 



.,^. 



k<l 



i-b-5 



£"^^,((0, 6)) = limlim / {G(x + ie)u,v) — {G(x — ie)u,v) dx 

' <5iO eio 27n j^^^ 

can come from the zeroes of the Wronskian. Put xi = 2{e + l)\/a'^ — 1, then for {a,b) containing 
precisely one of the xi 's we find 



Eu,v{{a,b)) = Eu,v{{xi}) 



1 

'2TTi 



M 



{G{z)u,v) dz. 



Note that the minus sign comes from the clockwise orientation of the rectangle with upper side 
X + ie, X £ [a, b], and lower side x — ie, x £ [a, b]. This residue can be calculated and we find 



Res 



1 



-Y^sin(7r(e — A)) sin(7r(— e — A)) 



Z = Xl 



.,^z 



2V~a 



Res 



Z=Xl 



sin(— evr + TTz/2\'a? — l) 



(-1) 



i+i 



vr 



■\/sin(7r(e — A)) sin(7r(— e — A)) 
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Since u'^{xi) = A^ {xi)vt{xi), because the zeroes of the Wronskian [u^ {z),v'^{z)] correspond to 
the zeroes of B^{z), see (4.4.8)| , we see that (pz and ^z are multiples of each other for z = xi. So 
we can symmetrise the sum in {G{z)u, v) again, and we find 

i-iy 



Eu,v{{xi}) 



nA {xi) 



\/sin(7r(e - A)) sin(7r(-e - X)){u, (j)^^ ){v,(j)x^). 



Using the explicit expression for A (xi) of (4.4.7) and the reflection identity for the T-function we 
arrive at 

'4(a2 - l)(a - V^^^^)\-efa + V^^^^\i 



Eu,v{{xi}) 



{u,(pxi){v,(t>xi) 



a + ^^2T7I ; \a- ^f^^^l' r(e + / - A)r(l + e + / + A) ' 

Note that in particular, (^^i are eigenvectors of L for the eigenvalue x/, and moreover these 
eigenspaces are 1-dimensional. Hence E{\xi\)(^xi ~ '^xr Recall Euy[{xi}) = {E(^{xi})u,v), and 



take u = V 



'JXi 



to find 



^Xl I 



4(a2 - l)(a - V^^^T)\ -^ /a + ^/^fl~^\l 



'Xl\ 



a + ^/^2-^T J Va _ ,/^r^iJ T{e + l- A)r(l + e + / + A) 

from which the squared norm follows. Using E[{xi})E{{xm}) = 5i^rriE{{xi}) and self-adjointness 
of E gives 

{4>xiAx^) = {E{{xi])(l)^i,E{{x„,])(t)xm) = {E{{Xm])E{{xi])(t)xiAxm) 
= 5l^m{E{{xi})4)xi,(t)xJ) = 5l^rn{(t)xiAxm) 

which proves the orthogonality. D 

(4.4.10) The orthogonality relations arising from Theorem (4.4.9) can be worked out and they 
give 

T{k + X + e + l)T{k + e-\) 



E 

fc=— oo 



(4(a2 - l)f 



T{k + l-l) 
1 



r(A; + l-p) 



2F1 
2 El 



k + e + \ + l,k + e-\ 1 



k + l-l '2 2^/^2rn: 

A; + e + A + l,A; + e-Al a 

k + l-p '2~27P^ 

4(a2 - l)(a - V^?^^)Yfa - V^^^T^ l 



^_^ ^_^ r(e + /-A)r(l + e + / + A). 

Since Theorem (4.4.9) also gives the completeness of this set of vectors, the dual orthogonality 
relations also hold, i.e. 



E 

/=— 00 



+ ^/^^^T\' 



1 



V^^^T/ T{e + l- A)r(l + e + / + A) 



1 



r(A: + l-/) 
1 



2 El 
2 El 



k + e + X + l,k + e- X 1 



k+l-l '2 2^/^2-^T 

m + e + A + l,m + e — A 1 a 

r(m + l-/)^^^V m + l-l 

(4(a^ - 1))' 



'2 2^/^^^n: 



Jk,m 



T{k + X + e + l)T{k + e - X) 



ija^ - l)ia - V^^^^)Y 
a + Va2 - 1 
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(4.4.11) Exercise. Use the results of [17| in order to calculate explicitly the spectral measures 
for the case \a\ < 1. In this case the spaces S^ are spanned by different solutions depending on the 
sign of Qz. 

In case |a| < 1 the situation changes considerably, and let us state briefly some of the results 
needed for the Meizner-Pollaczek case, i.e. \a\ < 1. The analogue of Lemma (4.4.4) is proved in 
entirely the same way, and we find that 



^„.,.,,....,)-.yi(^±J ^^r(^.e-., 



X 2F1 



X 2F1 



k + l + e^iz 'i_e±2iV'_ 

. Vr(l - A: + A - e)r{-k - e - X) 



^ ^ ' ' ^' ' "' T{l-k-e^iz) 

1 — k + X — e,—k — e — X 1 



l-k-e^iz '1- e=f 2iV' ^ 

satisfy the recurrence relation 

{2 sin ■ip)zukiz) = ak{cos7p,X,e)uk+i{z) + bk{cosip,X,e)uk{z) + ak^i{cosip,X,e)uk-i{z) 

with the explicit values for a^ and b^ of |(4.4.1)| with a = cosip, < ip < tt. This gives four linearly 
independent solutions of the recurrence relation. The connection formulae for these four solutions 
are given by 

U+{z) = A+{z) V+{z) + B+{z) V,-{z), U^{z) = A-{z) V+{z) + B-{z) V,~{z) 

and it follows easily from the explicit expressions and the assumptions on A and e of |(4.4.1)| that 
Uf{z) = U^{z) and V^{z) = V^{z), impl ying th at A-{z) = B+{z) and B-{z) = A+{z). The 
same formula for hyper geometric series as in K4.4.7)| can be used to find the connection coefficients; 

^ > ^ ^' T{iz-£)T{l+e-iz) 

B+(.) = (2..n^)''e'-0-i)^n_-^ " ^f' + ^ - ^J^d + ^ + .)r(. - A)^ 

r(A + 1 — iz)T[—X — iz) 

The asymptotic behaviour can be determined as in K4.4.5)| , and we find that (j)z = U^{z) for 
3=z > and 0^ = U-{z) for 3=z < and ^^ = V+{z) for 9z > and ^^ = V'iz) for 9z < 0. The 
explicit asymptotic behaviour can be used as in |(4.4.8)| to find the Wronskian 

[V-{z),V+{z)] = -i(2sinV)^"^"e-2^(^-5), 

and from this expression and the connection coefficients we can calculate all Wronskians needed. 
In order to find the spectral measure we have to investigate the limit e | of {G{x + i£)u, v) — 
{G{x — ie)u,v), X G M, and for this we consider, using the connection formulae, 

y+(x)c/+(x) Vj:{x)ur{x) 



B+{x)[V~{x), V+{x)] A-{x)[V+{x), V-{x)] 
V+{x)Vrix) + V,:ix)V+{x) A+{x)A-ix)Vi+ix)V+ix)+B-{x)B+ix)V^-ix)Vi-ix) 
[V-{x),V+{x)] ^ A-{x)B+{x)[V~{x),V+{x)] 
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which is obviously symmetric in k and /. Hence we can antisymmetrise the sum for the spectral 
measure and we find for u,v £ i'^{Z,) 

2TTi{u,v) = [(a-{x)B+{x){u,V+{x)){V+{x),v) + A-{x)B+{x){u,V^{x)){V'{x),v) 

+ A+{x)A-{x){u,V'{x)){V+{x),v) + B+{x)B"{x){u,V+{x)){V~{x),v) 

dx 

^ A-{x)B+{x)[V-{x),V+{x)] 

{u,U-{x)){U-{x),v) + {A-{x)B+(x) - A+{x)B- (x)){u,V- {x)){V- {x),v) 

dx 

^ A-{x)B+{x)[V-{x),V+{x)] 

describing the spectral measure, where we have used the relations between C/^ (z) and V^ (x) for the 
second equality. So we see that the spectrum of the corresponding operator is M, and by inserting 
the values for the Wronskian and the connection coefficients we see that the spectral measure is 
described by the following integral; 



{u,v) = 

{2smi;y^~^'e-^^^^~^\u,U-{x)){U-{x),v) 



\r{X + l-ix)T{-X-ix)\^ 



2tt 7Rr(-e - A)r(l + A - e)r(l + A + e)T{e - A) 
|r(A + l-ix)r(-A-ix 
\ix-e)T{l + e-ix)\'^ 



+ (l- ';;.^ " i"^".7 )(2sin^)-^+^-eM^-|)(n,y~(x))(y-(x),.)dx. 
\ \T{tx — e)T{l + e — ix]\'^ J 



We see that the spectral projection is on a two-dimensional space of generalised eigenvectors. Note 
that the general theory ensures the positivity of the measure in case u = v, and we can also check 



directly that, under the conditions on A and e as in (4.4.1) , the second term in the integrand is 
indeed positive. 

4.5. Example: the basic hypergeometric diflference equation. This example is based on 
Appendix A in [|13|| , which was greatly motivated by Kakehi |1^] and unpublished notes by Koorn- 
winder. The transform described in this section has been obtained from its quantum SU{1, 1) 
group theoretic interpretation, see |l^, ||l^ for references. On a formal level the result can be 
obtained as a limit case of the orthogonality of the Askey- Wilson polynomials, see [14| for a precise 



formulation. The limit transition descibed in (4.5.10)| is motivated from the fact that the Jacobi 



operators in this example and the previous example play the same role. 
(4.5.1) We take the coefficients as 



where we assume < q < 1, and r < 0, c > 0, d G M. This assumption is made in order to get the 
expression under the square root sign positive. There are more possible choices in order to achieve 
this, see |13, App. A]. Note that a^ and b^ are bounded for A; < 0, so that J" is self-adjoint. Hence, 



the deficiency indices of L are (0,0) or (1, 1) by Theorem (4.2.2) . 
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(4.5.2) In order to write down solutions of Lu = zu we need the basic hypergeometric series. 
Define 

fc— 1 n 

(a;g)fe = JJ(1 -ag'), k eZ>oU{oo}, (ai, ...,«„; g)A: = JJ(aj; g)fc, 

and the basic hypergeometric series 

fa,b \ ^ {a,b;q)k k 

The radius of convergence is 1, and there exists a unique analytic continuation to C\[l,oo). See 
Gasper and Rahman for all the necessary information on basic hypergeometric series. 

Lemma. Put 

,„2 _ ,2fc(c9^"Vrf^^;^)oo 



-•k 



(gi ^/r;g)c 



r^ t \ (^ \-k fdy,qdy/c q^-''c\ ^^ -N 



qy"^ ' ' (Pr 

then, with z = fi{y), we have that Uk{z) = Wkfkipiy)), Uk{z) = WkgkifJ-iv)), Uk{z) = WkFk{y) and 
Ukiz) = WkFk{y~^) define solutions to 

zuk{z) = akUk+i{z) +bkUk{z) + ak-iUk-i{z). 
Proof. Put Uk{z) = WkVk{z), then Vk{z) satisfies 

— h , 1 — i? 

2z Vk{z) = {d -— ) Vk+i{z) + q~''—: — Vk{z) + (d^^ -— ) Vk-i{z) 

ar dr dr 

and this is precisely the second order (/-difference equation that has the solutions given, see |^, 

exerc. 1.13]. D 



(4.5.3) The asymptotics of the solutions of Lemma (4.5.2) can be given as follows. First observe 
that W-k = d~^ as k —I- oo, and using 

2 k{r<i^^d^r/c,cq/d^r;q)oo _^ .^, i^^ , 

Wk = c ^ f^ ^ ^Wk = 0{C2 ), k^oo. 

{d^rq''/c,r,q/r;q)oo 

Now fkifJ-iv)) = 0(1) as fc ^ oo, and gkil^-iv)) = 0{{q/c)^) as k ^ oo. Similarly, F_k{y) = 
0{{dyf) ask^oo. 

Proposition. The operator L is essentially self-adjoint for < c < q^, and L has deficiency 
indices (1, 1) for q^ < c < 1. Moreover, for z G C\M the one- dimensional space S^ is spanned by 
wF{y) with ^{y) = z and \y\ < 1. For < c < q^ the one- dimensional space S^ is spanned by 
wf{z), and for g^ < c < 1 the two-dimensional space Sf is spanned by wf{z) and wg{z). 

Proof. In (4.5.1)1 we have already observed that the deficiency indices of L are (0, 0) or (1,1). Now 
2afc = q~^ y/c/dPr^ — ^(r + d'^r/c) + 0{q^), k —>■ oo, shows that the boundedness condition of 
Proposition |(3.4.6)| (ii) is satisfied if the coefficient of q~^ in a^ + Ok-i ± bk is non-positive. Since 
c > 0, dr G M, this is the case when (1 + q)\^ < c-\- q. For < c < q^ the inequality holds, so that 
by Proposition (3.4.6)K ii) also J+, and hence L by Theorem (4.2.2) , is essentially self-adjoint. 
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From the asymptotic behaviour we see that wf{z) and wg{z) are hnearly independent solutions 
of the recurrence in Lemma (4.5.3) , and moreover that they both belong to 5^ for q^ < c <\. The 
other statements follow easily from the asymptotics described above. D 



(4.5.4) The Wronskian 

[wF{y),wF{y-^)\ = lim akWk+iWk{Fk+i{y)Fk{y-^) - Fk{y)Fk+i{y^^)) = \{y-^ - y) 



using a^i ^ -^ as, k ^ —oo and the asymptotics of (4.5.3) . Note that the Wronskian is non-zero for 
y 7^ ±1 or z 7^ ±1. Since wF{y) and wF{y~^) are linearly independent solutions to the recurrence 
in Lemma 1(4.5.2) for z G C\M, we see that we can express fk{fJ-{y)) in terms of Fk{y) and Fk{y~^). 



These solutions are related by the expansion 

fkiKy)) = c{y)Fk{y) + c{y-^)Fk{y-^), 
, ^ ^ {c/dy, d/y, dry, q/dry; q)^^ 
(y-2,c,r,g/r;g)oo 

(4.3.2)]. This shows that we have 
[wfii^Liy)), wF{y)] = lc{y-'){y - y-'). 



(4.5) 



for c ^ Q' ^>o, y^ g^, see 



(4.5.5) Let us assume first that < c < g^, so that L is essentially self-adjoint. Then for z € 
we have 4>z = wf{z) and $2 = wF{y), where z = ii{y) and \y\ < 1. In particular, it follows that 
(px±ie -^ 4>x 8iS e I 0. For the asymptotic solution ^^ we have to be more careful in computing the 
limit of z to the real axis. For rr G M satisfying |x| > 1 we have ^x±ie ~^ '^^y as e J, 0, where 
y G (—1, 1)\{0} is such that /i(y) = x. If x G [—1, 1], then we put x = cosx = /"(e*-^) with x £ [0, vr]. 



and then (^x- 



wF„i^ and $ 



e'x 



x+ie 



wF^-ix as £ I 0. 



(4.5.6) We calculate the integrand in the Stieltjes-Perron inversion formula of (2.3.4)| using 
Proposition (4.3.3) and ( |4.3| ) in the case |x| < 1, where x = cosx = ^(e*"^)- Foi' u,v ^ 2?(Z) we 
have 



lim(G(x -|- ie)n, v) — {G{x — ie)u, v) 



lim > 



ej.0 



k<l 



i^x+ie)k{(px+ie)l i^x- 



^x+ie 



, ^x+ie] 



ie)k\SPx—ie)l 
^x—iet ^x—ie\ 



{uivk + UkVi){l - -6k,i) 



^E 



WkFkie '■^)wifi{cosx) WkFkie'^)wifi{cosx) 



k<l 



c(e*'^)(e-*x -e»x) 



c(e-*x)(e«x-e-»x) 



{uiVk + UkVi){l - -6k,i) 



2^\^WkWifi{cosx) 

k<l 

2E 



c(e^x)c(e-x)(e-x-e^x) (^'"^ + ^'="0(1 " ^^k,) 



k<l 



WkWifi{cosx)fk{cosx) / _ , _w. 1, ^ 
c(e*Xjc(e x)[e «x — e*xj ^ ' 2 



^ wifi{cosx)ui ^ Wkfk{cosx) 



c(e*x)c(e-*x)(e-»x - e«x) 



Vk 



using the expansion (^]^) and the Wronskian in (4.5.4). Now integrate over the interval (a, b) with 
— 1 < a < 6 < 1 and replacing x by cosx, so that ^dx = (e*^ — e~*^)dx/47r we obtain, with 
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a = COS Xa, b = COS Xb, and < X6 < Xa < vr, 

1 Z"-^" 
Eu,v{{a,b)) 



dx 



2^ 1^ (•^^)(cosx)(-^^)(cosx)|^^^,^^|2. 



Xb 



[J^u){x) = {u,(p^) = ^ wifi{cosx) 



Ui. 



l=—oo 



This shows that [—1, 1] is contained in the continuous spectrum of L. 



(4.5.7) For |x| > 1 we can calculate as in (4.5.6) the integrand in the Stieltjes- Perron inversion 
formula, but now we use that x = fj,{y) with \y\ < 1, see (4.5.5) . This gives 



lim(G(x + ie)u,v) 

eJ-O 



2E 

k<l 



WkFk{y)wifi{y) 
c{y~^){y - y~^) 



'uiVk + UkVi){l - -6k,i), 



and since the limit lime|o(G'(x + ie)u, v) gives the same result, we see, as in the case of the Meixner 
functions, that we can only have discrete mass points for |x| > 1 in the spectral measure at the 
zeroes of the Wronskian, i.e. at the zeroes of y i— > c{y^'^) with \y\ < 1 or at y = ±1. Let us assume 
that all zeroes of the c-function are simple, so that the spectral measure at these points can be 
easily calculated. 

The zeroes of the c-function can be read off from the expressions in ([4.5|), and they are {cq /d \ 
k G Z>o}, {dq'^ I k G Z>o} and {q^ /dr \ k G Z}. Assuming that \c/d\ < 1 and \d\ < 1, we see that 
the first two sets do not contribute. In the more general case we have that the product is less than 
1, since the product equals c < 1. We leave this extra case to the reader. The last set, labeled by 
Z always contributes to the spectral measure. Now for u,v & 2?(Z) we let Xp = fJ-iyp), yp = q^/dr, 
p G Z, with \qP/dr\ > 1, so that, cf. |(4.4.9)| , 

'{{xp}) = ReSj,=j;-i ( ^i^^_i^^] wkFk{yp^)wifi{xp){uiVk + UkVi){l - -5k,i) 



Em 



' y=yp \c{y- 

fj.{y)- Now from 



)y^ ■' '2 

after substituting x = n{y). Now from ( [4. 5] ) we find fk{xp) = c{yp)Fk{yp^), since c{yp^) = and 
we assume here that c{yp) ^ 0. Hence, we can symmetrise the sum again and find 

1 



Eu,v{{xp}) = Res, 



{Tu){xp)[j^v){xp) 



^-^-c(y-i)c(y)y. 
switching to the residue at yp. 

(4.5.8) We can combine the calculations in the following theorem. Note that most of the regularity 
conditions can be removed by continuity after calculating explicitly all the residues. The case of an 



extra set of finite mass points is left to the reader, cf. (4.5.7) , as well as the case of other choices 
of the parameters c, d and r for which the expression under the square root sign in o^ in (4.5.1) is 
positive. See |13, App. A] for details. 



Theorem. Assume r < 0, < c < g^, d G M with \d\ < 1 and \c/d\ < 1 such that the zeroes of 
y H^ c(y) are simple and c{y) = implies c{y~^) ^ 0. Then the spectral measure for the Jacobi 
operator on ^^(Z) defined by (4.5.1)| is given by, A cM. a Borel set, 

dx 



{E{A)u,v) 



{J='u) {cos x){^v) {cos x)- 



c(e^>^)p 



peZ,\qP /dr\>l,fj,{qP /dr)<=A 



cosxe[-l,l]nA 

Resy=qp/dr . _i. . . j{j'u){fi{qP/dr)){Tv){n{qP/dr)). 



Proof. It only remains to prove that ±1 is not contained in the point spectrum. These are precisely 
the points for which F{y) and F{y~^) are not linearly independent solutions. We have to show 
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that 

0, 



^±1 ^ (^)) ^-iid this can be done by determining its asymptotic behaviour as /c — > — cxd, see 
Ol for more information. D 



Take A = M and u = Ck and v 
2971-series as in Lemma (4.5.2); 



„ /fc(cosx)//(cosx)|^^^,^^| 



e;, then we find the fohowing orthogonaUty relations for the 
dx 



+ 



E ( 



1 



•'?^^.. .'7^^ ^k, 



R^^.=.-/^^ ,(^-i),(^)j A(M^))/K^(^)) = T^- 



tt't 



(4.5.9) In Theorem (4.5 



\qP/dr\>l '" ' ""'' "■ "■ ^ 

we have made the assumption < c < g^ in order to have L as an 



essentially self-adjoint operator. From the general considerations in (4.3.1) - (4.3.4) it follows that 
the previous calculations, and in particular Theorem (4.5.8) , remain valid for q^ < c < 1 if we can 
show that there exists a self-adjoint extension of L satisfying the assumptions of (4.3.2). It suffices 
to check (2) of (4.3.2), or, by Lemma |(4.2.3) , the existence of a 6 £ [0, 2-k) such that 



lim [wf{z),e'^wF{i{l - V2)) + e-'^wF{-i{l - V2))]k = 0. 



(4.6) 



Indeed, wF{±i{l — \/2)) is the element ^±i up to the normalisation of the length. This is not 
important for showing the existence of 9 satisfying ( |4.6| ). 

For (|4.6| ) we need to know the asymptotic behaviour of Fk{y) as A; ^ 00. The same result in 
basic hypergeometric series that results in (5^) can be used to prove that 



a{y) 



Fk{y) = a{y)fk{y) + Ky)9k{y), 

{qdy/c, qy/d, qcy/dr, dr/cy; q)oo 



b{y) 



{dy, cy/d, q^y/dr, dr/yq; q)^ 



{qy'^,q/c,qc/d'^r,d'^r/c;q)oo ' {qy^,c/q,qc/d'^r,d'^r/c;q)c 

for y^, c ^ q^, so that Ffc(y) = a{y)0{l) + b{y)Oi{q/c)''). It follows that 

1, 



lim [wf{z),wF{y)]k , , 

fc^oo 2 dr 



1 



q 



'-)b{y), 



so that the limit in ( [4.6[ ) equals 

l:fl(l- 

\ dr \^ 



■ {e'%{i{l - V2)) + e-'%{-i{l -V2))). 



2\\Fiiil - V2))\\ 

The term in parentheses is 2K(e 6(i(l — v2))), which is zero for 9 = — ^ + arg6(i(l — v2))- 
(4.5.10) The Meixner functions can formally be obtained as a limit case as g | 1 from the spectral 
analysis of the second order g-difference operator as considered here. For this we make the following 
specialisation; c 1-^ gs"^, d 1-^ g^"*"'^s~^, r h^ g~^^^, where e and A have the same meaning as in 
(4.4.1) . Note that r < is no longer valid but the operator is well-defined under suitable conditions 



on A, cf. (4.4.1) . The operator L with (4.5.1) has the same type of spectral measure (exercise). 
Now consider the operator 



L„ 



2L 



1 



Lqfk 



4 /fc+1 + blfk + aLi fk-i 



for fk = B-k, then a calculation shows that a^ ^ a/c, ^^ — > bk with a^ and b^ as in (4.4.1) with 
2a = s + s^^ as g I 1. We now assume a > 1 as in the previous subsection on Meixner functions. 

Now the operator Lg has continuous spectrum supported on [(—2 — s — s~^)/(l — q),{2 — s — 
s~^)/(l — q)]. For s + s~^ > 2 the continuous spectrum will disappear to —00 as g j 1; it will tend 
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to (—00, 0] for s + s~^ = 2, and it will tend to M if < s + s~^ < 2. The discrete spectrum (or at 
least the infinite number of discrete mass points) is of the form 

p-l+e _ 1 n^~'^~-P — 1 

s^ + s-^^—^ -, peZ, \sqP+'''\ > 1. 

1-q 1-q 

As q 1 1 this tends to {p + e — l){s~^ — s) with p G Z for |s| > 1 and it will disappear for |s| < 1. 



So for a > 1 only the discrete spectrum survives and this corresponds precisely to Theorem (4.4.9 
For < a < 1 the discrete spectrum disappears in the limit. 

It is also possible to show that the solutions for the g-hypergeometric difference operator tend 
to to solutions of the operator for the Meixner functions. This requires the transformation as in 



(4.4.6), and one of Heine's transformations. 

So we have motivated, at least formally, that the orthogonality relations for the Meixner functions 
can be obtained as a limit case of the 2(/'i-series that arise as solutions of the second order g-difference 
hypergeometric equation, and that the limit transition remains valid on the level of (the support 
of) the spectral measure. For the Laguerre and Meixner-Pollaczek functions the limit remains valid 
on the level of the support of the spectral measure. 
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